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ABSTRACT 

The  general  topic  of  the  research  reported  in  this  thesis  is  the 
study  of  nonunifonn  block  sampling  theory  and  its  application  to  the 
study  of  correlation  functions  of  sampled  signals.  The  research 
evolved  from  the  specific  problem  of  trying  to  reduce  the  odd  harmonic 
terms  which  appear  in  the  correlation  function  of  hard  limited  and 
uniformly  sampled  sinusoidal  signals.  These  odd  harmonic  terms  arise 
because  of  the  synchronism  between  the  periodicity  of  the  uniform 
samples  and  the  periodicity  of  the  odd  signal  harmonics  produced  by 
hard  limiting. 

The  general  approach  to  the  problem  is  to  use  nonuni  fOrm  sampling 
to  reduce  the  undesirable  periodicity  effects  and  thus  to  reduce  the 
odd  harmonic  terns  in  the  correlation  function.  A  particularly  useful 
class  of  nonuni  form  sampling  functions  for  this  application  is  shown  to 
be  “block  sampling  functions;"  l.e. ,  short  sequences  (or  blocks)  of 
nonuniftorn  samples  which  are  repeated  periodically  in  time. 

As  a  preliminary  to  a  study  of  block  sampling  functions  and 
their  application  to  the  odd  harmonic  problem,  a  detailed  study  ox  the 
general  properties  of  correlation  functions  for  large  TV  signals  with 
uniform  tangling  is  made.  These  results  are  then  <u tended  to  the  case 
of  sinusoidal  signals  with  uniform  stapling  to  illustrate  the  role  of 
both  clipping  and  sampling  in  the  generation  of  the  undesirable  odd 
harmonic  tense  in  the  correlation  function.  The  uniform  sampling 
correlation  function  results  along  with  the  results  of  the  block 
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sampling  spectrum  analysis  then  fora  the  basis  for  analysis  of 
correlation  functions  of  sinusoidal  signals  with  nonuni  fora  block 
sampling. 

In  the  study  of  correlation  functions  of  large  TV  signals,  two 
alternative  definitions  of  sampled  correlation  functions  are  compared: 
the  regular  correlation  function,  and  a  quadrature  correlation  function 
consisting  of  an  average  of  a  regular  correlation  function  and  the 
correlation  function  of  the  Hilbert  transforms  of  the  correlated 
signals.  By  eliminating  the  second  harmonic  sampling  effects,  the 
quadrature  correlation  function  is  shown  to  achieve  sore  accurate 
correlation  for  deterministic  signals  and  lower  correlation  variance 
for  random  signals.  A  comparison  of  the  saapled  quadrature  cc. 1 elation 
function  variance  for  clipped  and  undipped  narrowband  random  processes 
shows  that  the  variance  for  the  clipped  signals  could  be  reduced  by  as 
much  as  3  dB  over  the  undipped  signal  varlanca  by  sampling  at  a 
sufficiently  high  rate. 

In  extending  the  study  of  quadrature  correlation  Auction 
properties  to  sinusoidal  signals,  uniform  sampling  la  shown  to  sake 
tbs  correlation  Auction  periodic  (as  s  Auction  of  frequency)  with  s 
repetition  period  equal  to  the  sampling  rate,  while  clipping  is  traduce  s 
extraneous  correlations  for  input  signal  frequencies  which  are  odd 
multiples  of  the  reference  signal  frequency.  Vhee  both  clipping  sod 
ualfern  sampling  are  employed,  each  odd  harmonic  correlatloe  Auction 
reepooee  is  repeated  at  period  frequeecy  Intervals  which  a-e  related 
to  the  aanplU*  rate,  thus  giving  rise  te  the  odd  harmonic  problem  for 
ear  re  let  ion  Auction*  of  eimeaoidnl  signals. 
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After  proposing  nonuniform  block  sampling  as  &  possible  solution 
to  the  odd  harmonic  problem,  the  frequency  spectrum  properties  of 
ideal  block  sampling  are  studies  in  detail.  These  properties  include 
amplitude  and  phase  characteristics  of  each  spectral  component, 
spectral  periodicity  conditions,  conservation  of  spectral  energy 
properties,  and  conditions  for  minimum  variance  of  the  sampling  power 
spectrum.  As  an  example  of  sampling  spectrum  amplitude  shaping  by 
choice  of  sample  times,  a  block  sampling  function  example  is  presented 
which  uses  a  pseudo-random  number  generator  for  sampling  time  selection. 
This  sampling  technique  achieves  a  unique  spectrum  shaping  by 
suppressing  the  spectral  amplitude  for  small  values  of  frequency,  a 
property  that  is  shown  to  be  important  in  the  reduction  of  odd 
harmonic  responses  in  correlation  functions  of  clipped  sinusoidal 
signals. 

Considering  the  quadrature  correlation  function  of  two  clipped 
sinusoidal  signals  v.th  block  sampling  as  a  function  of  the  two 
sinusoidal  frequencies,  it  is  shown  that,  although  tha  odd  harmoni- 
effects  arc  greatly  reduced,  certain  "Vorst  Case"  combinations  of  the 
two  sinusoidal  frequencies  result  la  correlation  functions  whose 
average  value  is  greater  than  those  fUr  all  other  frequency  conbiaa- 
t ions,  la  addition  to  deriving  the  critical  frequency  relationships, 
sa  analytic  expression  Is  attained  relating  tbs  correlation  function 
far  these  Verst  Case  frequencies  to  the  cos  me  lasts  of  tbs  block  sampling 
frequency  apectnse.  teas  experimental  renal ts,  based  no  s  digital  eeqwur 
slmslatlaa,  are  shown  Vs  verify  the  analytical  aadel.  It  is  concluded 
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that  by  proper  choice  of  sampling  parameters,  block  sampling  can  reduce 
the  undesirable  odd  harmonic  correlation  function  responses  due  tc 
*  »  clipping  and  sampling  by  any  desired  degree. 
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CHAPTER  I 


INTRODUCTION 

1.1  General  Statement  of  the  Problem 

The  general  topic  of  tne  research  reported  in  this  thesis  i6  the 
application  of  nonur.iform  sampling  theory  to  the  study  of  correlation 
functions  of  sampled  signals  that  are  limited  to  a  band  of  frequencies. 
In  addition  to  nonuniform  sampling,  the  research  considers  effects  of 
uniform  sampling  and  hard  limiting  on  the  correlation  functions  of 
signals  with  large  TW  products  as  well  as  sinusoids.  The  research 
evolved  from  the  specific  problem  of  trying  to  find  a  method  to  reduce 
the  odd  harmonic  terms  which  appear  in  the  correlation  function  of  a 
sinusoidal  signal  which  has  been  hard  limited  and  uniformly  sampled. 
These  odd  harmonic  terms  arise  because  of  the  synchronism  between  the 
periodicity  of  the  unifo  m  samples  and  the  periodicity  of  the  odd 
signal  harmonics  produced  by  hard  limiting.  This  research  considers 
the  use  of  a  special  class  of  nonuniform  sampling- -"block  sampling;' 
i.e.,  nonuniform  samples  with  a  periodically  recurrent  pattern--to 
eliminate  this  synchronism  and  reduce  the  odd  harmonic  correlation 
function  terms. 

1.2  Importance  of  the  Study 

The  correlation  function  is  an  important  analysis  tool  in  a 
large  number  of  applications  £5,  26,  31,  32].  It  is  particularly 
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useful,  for  example,  in  the  study  of  the  response  of  linear  systems  to 
random  inputs  [33]  where  deterministic  analysis  methods  do  not  apply 
and  in  determining  the  statistics  of  signals  in  the  presence  of  noise 
[13].  Since  the  correlation  function  and  the  power  spectrum  are 
Fourier  transforms  of  each  other,  it  is  also  important  in  spectral 
analysis  [13].  Other  applications  include  Hilbert  transform  theory 
[llj]  and  the  correlation  function  method  [H]  for  finding  roise 
statistics  at  the  output  of  certain  no.ilinear  devices. 

The  operations  of  sampling  and  a  hard  limiting  (clipping)  have 
become  conmon  in  many  physical  systems  [29*  26,  lit ,  35].  These 
operations  are  often  necessary  because  of  physical  limitations  which 
make  available  only  sampled  information  o*-  signal  polarity  information. 
Sampling  has  become  an  important  item  in  modern  communication  systems 
since  the  advent  and  wide  spread  use  of  digital  techniques.  Examples 
of  the  use  of  sampling  include  pulse  code  modulation  (FCM)  systems  and 
pulsed  radar.  In  most  sampling  applications  in  the  literature,  the 
concern  is  reconstruction  of  the  original  signal  from  its  samples  [l!»]. 
The  basis  for  representing  continuous  signals  by  sampled  values  taken 
at  discrete  time  instants  lies  in  the  well-known  sampling  theorems  [It], 
For  signals  limited  to  a  frequency  band,  the  sampling  theorems  are 
usually  presented  in  terms  of  quadrature  sampling  [20];  i.e.,  taking 
samples  of  both  the  signal  and  its  Hilbert  transform  at  a  rate 
determined  by  the  signal  bandwidth.  In  this  thesis,  we  are  not 
concerned  with  the  reconstruction  of  the  original  signal  from  its 
sampled  values,  but  in  the  effects  of  sampling  on  the  auto-correlation 


3 


function  of  a  signal  or  on  its  cross-correlation  function  with  other 
signals . 

A  companion  process  to  sampling  in  many  digital  applications  is 
that  of  signal  quantizing.  Many  communication  systems  require  that 
signal  samples  be  converted  to  digital  numbers  at  some  point  within 
the  system  so  they  can  be  manipulated  by  digital  computers  or  other 
digital  devices.  An  important  special  case  of  quantizing  is  hard 
limiting  [37,  23],  or  clipping,  in  which  only  two  levels  of 
quantization;  i.e.,  signal  polarities  are  employed.  In  this  research, 
we  will  consider  the  implications  of  hard  limiting  on  the  correlation 
function  of  sampled  signals  that  are  limited  to  a  beuid  of  frequencies 
[30,  38], 

Although  most  of  the  sampling  theory  work  in  the  literature  is 
concerned  with  uniform  sampling,  many  important  applications  [36,  28] 
either  intentionally  [1]  or  unintentionally  [2,  9]  madce  ure  of 
nor.uniform  sampling.  For  this  reason,  a  good  deal  of  research  [ill , 

2h ,  5,  21,  J»0]  has  been  done  on  the  problem  of  signal  reconstruction 
from  nonuniformly  spaced  samples  of  a  signal.  The  effect  of  nonunifcrm 
sampling  on  correlation  function  properties  of  sinusoidal  signals  is 
to  be  a  major  consideration  in  this  thesis. 

1.3  Previous  Related  Studies 

The  special  class  of  nonuniform  sampling  functions — those  with 
a  periodically  recurrent  pattern- -applied  In  this  thesis  to  the  problem 
of  reduction  of  odd  harmonics  in  correlation  functions  of  clipped 
sinusoidal  signals  has  been  considered  previously  by  several  authors 
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[4l,  2h,  36,  7].  Yen  [Ul],  using  the  term  "Recurrent  Nonuniform 
Sampling,"  and  Kohlenberg  [2k],  using  the  term  "p'h  Order  Sampling," 
both  consider  the  problem  of  reconstruction  (interpolation)  of  a 
continuous  signal  from  its  ncnuniformly  sampled  values.  Tou  [36], 
using  modified  z-transform  theory,  analyzes  the  response  of  sampled- 
data  systems  using  "Cyclic  Variable-Rate  Sampling,"  while  Beutler  and 
Leneman  [7]  consider  'Periodically  Recurrent"  nonuniform  sampling 
functions  as  a  special  class  of  processes  in  the  general  theory  of 
stationary  point  processes . 

More  specifically,  Yen  [U]  derived  a  reconstruction  formula 
for  low  pass  signals  (bandlimited  to  W  Hz)  which  were  sampled  by  a 
recurrent  nonuniform  sampling  function  with  N  samples  in  each 
periodic  sample  block  .  Kohlenberg  [2^]  derives  a  more  general  method 
for  obtaining  reconstruction  formulas  for  the  same  sampling  function, 
requiring  only  thaT  the  sampled  signal  possess  a  Fourier  Transform. 

While  Kohlenberg  s  method  applies  in  principle  to  signals  for  an 
arbitrary  number  N  of  samples  per  recurrent  sample  block,  he  applies 
it  only  for  N  =  2  obtaining  a  second-order  sampling  reconstruction 
formula.  It  is  apparent  that  his  method  becomes  very  cumbersome  for 
N  >  2.  Kohlenberg  expresses  recurrent  nonuniform  sampling  with  N 
samples  per  block  as  the  sum  of  N  uniform  sampling  functions  with 
different  time  origins  wnere  the  uniform  sampling  period  of  each  function 
is  equal  to  .he  block  interval  length  in  seconds.  This  representation 
proved  useful  in  the  study  of  Block  Sampling  in  this  thesis  research. 

Tou  [36]  models  nonuniform  block  sampling  in  a  slightly 
different  way,  shifting  the  continuous  signal  (in  parallel)  with  N 
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different  time  advances  (corresponding  to  the  time  of  each  sample 
instant  in  the  block),  sampling  all  N  signals  at  the  same  instant, 
delaying  each  sampled  signal  by  an  amount  equal  to  the  time  advance 
before  sampling,  and  finally  sunning  the  N  sampled  signals  to  form 
the  resultant  sampled  signal.  A  modified  version  of  Tou's  model  was 
found  useful  in  Chapter  V  for  computing  the  correlation  function  of  a 
clipped  sinusoidal  signal  with  nonuniform  block  sampling. 

Block  sampling  is  shown  by  Beutler  and  Leneman  [7]  to  be  a 
special  example  of  a  S.P.P  (stationary  point  process)  if  the  initial 
sampling  instant  is  allowed  to  be  u  random  variable.  Studying  the 
problem  of  analyzing  stochastic  sampling  of  wlde-sense  stationary 
random  processes,  Beutler,  Leneman,  and  Lewis  [6,  7»  27*  23]  formulate 
a  basic  theory  for  S.P.P.  deriving  statistics  on  the  number  of  points 
in  intervals  and  on  forward  and  backward  recurrence  times,  and  apply 
the  theory  to  a  number  of  problems  including  spectral  analysis  of 
randomly  modulated  processes  and  mean  square  error  in  t  te 
reconstruction  of  signals  from  randomly  timed  samples.  In  this  thesis, 
block  sampling  will  be  treated  as  a  deterministic  process,  making  use 
of  Fourier  transform  theory  to  derive  and  study  its  complex  spectrum 
properties. 

1.1»  Scope  and  Limitations  of  the  Study 

The  basic  topic  of  this  research — the  study  of  nonuniform  block 
sampling  and  itr  application  in  reducing  the  odd  harmonic  terms  in 
correlation  functions  of  clipped  sinusoidal.  signals--is  covered  in 
Chapters  I”  and  V.  In  Chapter  IV,  nonuniform  '"block  sampling"'  is 
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defined  and  its  spectral  properties  are  derived.  These  results  are 
then  applied  in  Chapter  V  to  the  odd  harmonic  problem  in  correlation 
functions.  Experimental  results  from  a  digital  computer  simulation 
are  include .  to  verify  some  of  the  theoretical  results . 

As  a  preliminary  to  the  nonuniform  sampling  work  of  Chapters  IV 
and  V,  a  study  of  the  general  properties  of  correlation  functions  for 
large  TV  signals  including  the  effects  of  uniform  sampling  and 
clipping  is  made  in  Chapter  II.  These  results  are  then  extended  in 
Chapter  III  to  quadrature  correlation  functions  of  sinusoidal 
signals  illustrating  the  role  of  both  clipping  and  uniform  sampling 
in  the  generation  of  the  odd  narmor.^c  correlation  function  terms. 


CHAPTER  II 


CORRELATION  FUNCTION  PROPERTIES  OF  SIGNALS  WITH  LARGE  TW 
THAT  ARE  LIMITED  TO  A  BAND  OF  FREQUENCIES 

2.1  Definition  of  Regular  and  Quadrature  Correlation  Functions 
In  this  chapter,  a  study  of  the  general  properties  of 
correlation  functions  of  signals  with  large  time-bandwidth  (TW) 
products  will  be  made.  After  comparing  two  alternative  definitions 
for  correlation  functions  of  uniformly  sampled  signals,  the  effects 
of  clipping  will  be  considered. 

The  two  alternative  correlation  function  definitions  to  be 
compared  are  the  regular,  or  standard,  sampled  correlation  function 
(20)  and  the  quadrature  correlation  function  which  makes  use  of 
Hilbert  transform  pairs  of  signal  samples. 

The  regular  correlation  functions  * (MT  )  and  ^(MT  )  are 

S  8 

defined  by 

M 

*<”.>  -  R  ?  x(mTfl)  y(mTg)  (2.1) 

m-1 

and 

M 

»(MTg)  =  i  x(mTg)  y(mTg),  (2.2) 

m=l 

where  T  is  a  uniform  sampling  period,  M  is  the  number  of  samples 
s 
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in  a  signal  duration  T,  x(t)  and  y(t)  are  real-valued  signals  which 
can  be  either  deterministic  or  random,  and  y(t)  is  the  Hilbert 
transform  of  y(t).  The  usual  correlation  function  dependence  on 
time  shift  t  can  be  found  by  setting  x(t)  in  Equations  (2.1)  and 
(2.2)  equal  the  time  shifted  version  of  itself;  i.e.,  by  setting 
x(t)  =  x(t  -  :).  The  envelope  of  the  regular  correlation  function 
can  be  written  as 

£.  (MTg )  «-W»2(MTs)  +  *2(MTs)  .  (2.3) 

An  analytical  model  for  the  regular  correlation  function  is 
shown  in  Figure  2.1  using  complex  signal  representation  [15].  In  this 
model,  the  sampled  pre-envelope  p  (mT  )  of  the  reference  signal 

y  s 

y(t),  defined  by 

py(MTs)  -  y(»Ts)  +  j  y(mTs)  , 

is  correlated  with  the  sampled  input  signal  x(mT  )  giving  the 

s 

following  complex  correlation  function: 

M 

*  H  Py("T.) 

m=l 

M  M 

*  k  x(mTs)  y(mTs)  +  *  k  x^mTs)  y(mT8) 

»-l  m=l 

«  •(MT  )  +  j  *(MT  ). 

O  O 


Flfur*  Z  1  AMljrtlc  Modal  of  >*gul*r  Correlation  Function 
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Thus,  the  complex  correlation  function  z(MT  )  contains  both 

s 

quadrature  components,  ®(MT  )  and  o(MT  ),  and  the  envelope, 

s  s 

&MTJ;  i.e., 

Re[z(MT  )]  =  $(MT  )  =  'in-phase1'  component  of  the  regular 

s  s 

correlation  function, 

Im[z(MT  )]  =  $(MT  )  -  quadrature  phase''  component, 

s  s 

and 

z(MT  )  -  ^(MT  )  -  regular  correlation  function  envelope. 

II  5  |  S 

We  will  next  consider  a  modified  version  of  the  regular 

correlation  function  model  of  Figure  2.1.  The  "Quadrature  Correlation 

Function  '  model  of  Figure  2  2  differs  from  the  Regular  Correlation 

Function  model  of  Figure  2.1  in  that,  for  the  quadrature  version, 

Hilbert  transforms  are  taken  of  both  the  input  and  reference  signals 

rather  than  just  the  reference  signal,  Instead  of  sampling  the  signal 

x(t)  every  T  seconds  as  shown  in  Figure  2  1,  the  quadrature 
s 

correlation  requires  pairs  of  samples,  x(2mT  )  and  x(2mT  ),  every 

s  s  * 

2Tg  seconds,  giving  the  same  total  number  M  cf  samples  during  the 

finite  signal  duration,  where  M  is  assumed  to  be  an  even  integer. 

The  resulting  equations  for  the  in-phase  component  ®.(MT  ),  tne 

s 

A  r* 

quadrature-phase  component  ®q(Mts)»  and  the  envelope  c^(MTs)  of 
the  quadrature  correlation  function  are 
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NOTE :  H  INDICATES  THE  OPERATION  OF  TAKING  THE  HILBERT  TRANSFORM. 


SAMPLER 


Figure  2.2  Analytic  Model  of  Quadrature  Correlation  Function 
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[x(2mTs)  y(2mTg)  +  x(2mTg)  y(2mTg)], 

(2.M 

[x(2mTs)  y(2mTg)  -  x(2mTg)  y(2mTg)], 

(2-5) 

£q(KTs)  *V*q2  (MTg)  +  *Q2  (MTg)«  (2.6) 

A  comparison  of  Equations  (2-M  and  (2  1)  show  that  both  types  of 
correlation  function  definitions  require  the  same  number  of 
numerical  computations  (products  and  sums)  per  second  to  compute  one 
value  of  the  correlation  function.  However,  successive  values  of  the 

quadrature  correlation  function  can  only  be  computed  at  2T  second 

s 

intervals  instead  of  at  Tg  second  intervals  if  the  time  delay  i 
is  performed  after  sampling. 

In  the  analytic  quadrature  correlation  function  model  shown  in 
Figure  2.2,  the  pre-envelope  of  the  reference  signal  y(t)  is 
correlated  with  the  complex  conjugate  of  the  pre-envelope  of  the  input 
signal  x(t)  giving  the  following  equation  for  the  complex  quadrature 
correlation  function,  zQ(MT  ): 

Vi  S 

M/Z 

W  •  s  <a»V  Py(2»Ts) 

m-1 


M/2 

i 

M  / 
m=l 


“  M 


nrl 
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M/2 

=  k 


M  (_ 

m=l 


[x(2mT  )  -  jx(2mT  ) ] [y(2raT  )  +  jy(2mT  )] 

5  5S  S 


M/2 


=  |  )  [x(2nTs)  y(2mTs)  +  x(2mTs)  y(2mTs)] 


nr=l 


4  J  i  '  [x(2mT.)  y(2mTs)  -  x(2mT#)  y(2mTg)] 


m=l 


*  *Q  (MTS)  ♦  J  \  (MT#). 


Taking  the  real  and  imaginary  parts  and  the  absolute  value  of  *Q(MT  ) 

Vi  S 

gives : 


Re  [*Q(MTg)l 
t*Q(MTs)] 

K(mV| 


£s(MTt), 


the  in-phase  component, 

the  quadrature-phase  component, 

the  correlation  envelope. 


2.2  Comparison  of  Regular  and  Quadrature  Sampled  Correlation 
Functions 

A.  Elimination  of  the  Second  Harmonic  by  ^adrature  Correlation. 
Consider  the  equations  for  the  two  basic  correlation  techniques  being 
compared,  regular  correlation  and  quad  *ature  correlation.  From 
Equation  (2.1),  the  regular  correlation  function  ’in-phase”  component 
is: 


ll* 


M 

•(“,)  •  s  :  *<»T,>  y("Tt) 

m=l 


(2.7) 


and,  from  Equation  (2.1* ),  the  quadrature  correlation  function 
"in-phase’  component  is: 


♦q(mts) 


1 

M 


M/2 


[x(2mT#)  y(2mTs)  +  x(2mTg)  J(2mTg)]. 


m^l 


(2.8) 


Equauions  (2.7)  and  (2.8)  are  both  sampled  approximations  to  the 
desired  analog  correlation  function  given  by 


•  (T) 


x(t)  y(t )  dt, 


(2.9) 


wnere  x(t)  and  y(t)  are  continuous  real-valued  functions  that  are 
approximately  limited  in  time  T  and  bandwidth  W  such  that 
TW  »  1  The  basic  question  that  we  now  wish  to  examine  is  how 
accurately  the  sampled  correlation  functions  of  Equations  (2.7)  and 
(2  8)  approximate  the  desired  continuous  correlation  Tunction  of 
Equation(2.9)  for  two  specific  cases  The  two  cases  considered  here 
are:  (l)  x(t)  and  y(t)  are  both  deterministic  functions  for 
which  Fourier  transforms  can  be  written  [as  a  special  case,  y(t) 
could  be  a  delayed  version  of  x(t)]f  and  (2)  both  x(t)  and 
y(t)  (either  deterministic  or  random)  can  be  represented  by 
narrow-band  representations  of  the  form  [37] 
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s(t)  =  a(t)  cos  a)Qt  -  b(t)  sin  a^t, 

where  a(t)  and  b(t)  are  low-pass  functions  and  cdo  is  the  center 
frequency  of  the  signal  s(t).  Case  1  thus  considers  the  sampling 
errors  in  obtaining  the  correlation  function  between  two  deterministic 
signals,  and  Case  2  considers  the  effects  of  sampling  on  the 
correlation  function  variance  when  a  reference  signa^  is  correlated 
with  a  narrow  band  random  process.  In  both  cases,  it  will  be  found 
that  the  sampled  quadrature  correlation  function  represents  a  more 
accurate  approximation  to  Equation  (2.9)  if  the  sampling  rate  is 
sufficiently  high  and  that  the  sampling  errors  associated  with  the 
regular  correlation  function  are  a  function  of  a  signal  bandwidth, 
center  frequency,  and  sampling  rate. 

The  basic  phenomenon  that  makes  quadrature  correlation  more 
accurate  than  regular  correlation  in  both  cases  is  that  quadrature 
correlation  eliminates  the  second  harmonic  in  the  frequency  spectrum 
of  the  unsampled  product  of  x(t)  and  y(t).  This  will  be  shown  by 
obtaining  expressions  for  the  Fourier  transforms  of  the  products 
[x(t )  •  y(t)]  for  regular  correlation  and  ^  [x(t)  y(t)  -1-  x(t)  y(t)j 
for  quadrature  correlation. 

Considering  Case  1,  let  X(u>)  and  Y(u))  represent  the  Fourier 
transforms  of  x(t)  and  y(t),  and  let  $(oo)  and  Y(to)  represent 
the  Fourier  transforms  of  x(t)  and  y(t).  Assume  x(t)  and  y(t) 
are  deterministic  signals  whose  non-negligible  positive  spectrums,  as 
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illustrated  in  Figure  2.3,  are  confined  to  the  frequency  interval 
(co  -  ^-)  <  co  <  (<jd  +4),  where  co  is  the  center  of  the  frequency 

O  L  O  b  O 

band,  A  -  2rtW,  and  W  is  the  signal  bandwidth  in  Hertz.  According 
a  A 

to  [15] ,  X(u))  and  Y(co)  can  be  expressed  in  terms  of  X(co)  and 

Y(co)  as 


X(co) 


I  -jXfco)  for  u)  >  0 
|  jX(co)  for  co  <  0  , 


(2,10) 


and 


a  I  -jY(oo)  for  co  >  0 

Y'oo)  = 

jY(co)  for  co  <  0  . 


(2.11) 


The  Fourier  transform  Fr(co)  of  the  regular  correlation 
product  x(t)  y(t)  is: 


FrM  =J[x(t)  y(t)]  =  x(t)  y(t)  e"^  dt.  (2.12) 


Substituting  the  expressions  for  the  inverse  transforms  of  x(t) 
and  y(t)  into  Equation  (2.12)  gives  FR(co)  as  the  convoluticn  of 
X(co)  and  Y(co): 


FR(u>) 


X(B)  Y(u>B)  d&. 

“*  —  oc 


(2.13) 


Likewise,  for  the  product  x(t)  y(t),  if 


F2(co)  =  !^[x(t)  y(t )]  =  J  x(t)  y(t)  e~J  ^  dt, 


then 
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i 


|X(„|| 


-2«»  -«p  0  2«e 


-2w0  -we  0  m0  2*10 


Figure  2.3  Typical  Non-Zero  Signal  Spectrum  Ranges 


r-l|CM 
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f2(u>)  =  X(e)  Y(u>-3)  dp.  (2.1U) 

-at 

The  spectrum  of  the  quadrature  correlation  product 

[x(t)  y(t)  +  x(t)  y(t)]  in  terms  of  FR(u))  and  F2(u>)  *s: 

Fft(u>)  =  [x(t)  y(t)  +  $(t)  y(t)]j- 

=  |  LFR(u>)  +  F2(u,)).  (2.15) 

An  examination  of  the  frequencies  at  which  FR(u))  and 
Fq(cu)  are  non-negligible  will  show  that,  for  F^(u)),  the  second 
harmonic  terms  around  u>  =  2u>  cancel  when  the  sum  of  F_,(a>)  and 

O  H 

F^ (u>)  is  taken  in  Equation  (2.15).  Consider  first  the  non-negligible 

regions  of  the  regular  correlation  spectrum  Fp(co).  Hie  graphical 

convolution  shown  in  Figure  2.U  of  X(uo)  and  Y(o>)  from  Figure  2.3 

shows  that  the  non-negligible  portions  of  F  (u»)  ere  given  by 

n 

-A  <  <d  <  A,  (modulation  term) 

(2uqo  -  A) <  u>  <  (2uq  +  A),  (positive  second 

harmonic ) 

(-2llo  -  A)  <  a)  <  ( “2u-0  +  A),  (negative  second 

harmonic ) . 

(2.16) 

The  corresponding  expressions  for  Fp(x)  obtained  from  Equation  (2.13) 


are: 
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-to  +  A/2 
i  i° 

Fp(u>)  =  i  /  X(p)  Y(co-a)  dp 

R  2n  J  -to  -  A/2 


x(p)  y(u-p)  da, 

for  -A  <  cu  <  A, 

x(a)  Y(o>-a)  da, 

for  (2u>o  -  A)  <  to  <  (2a>o  +  A), 

x(e)  Y(to-a), 

for  (-2coq  -  A)  <  to  <  (-2uo  +  A). 

(2.17) 

An  examination  of  F2(u.)  =  ^[x(t)  y(t)]  from  Equation  (2.1b) 

shows  the  same  non-negligible  spectral  regions  given  by 
Equation  (2.16),  where  F2(to)  can  be  derived  using  the  Hilbert 
transform  spectral  properties  of  Equations  (2.10)  and  (2.11)  as: 


♦i 


J 


to  +  A/2 


co  -  A/2 
o  ' 


1 

2n 


io  +  A/2 
JC0Q  -  A/2 


1 

2n 


-co  +  A/2 

J  , 

^  -co  -  A/2 
o  ' 


f2M 


i 

2x 


J 


-co  +  A/2 


-co  -  A/2 
o  ' 


to  +  A/2 

r 

JtOQ  -  A/2 


[jx(a)j  [-jy(u»-p)j  da 

(-jx(a)]  [jy(p)j. 


for  -A  <  to  <  A, 
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1 

2n 


+  A/2 
-  A/2 


[-jx(e)]  [-jY(o>p)i  dp, 


for  (2a>o  -  A)  <  co  <  (2coq  +  A), 


1 

2n 


-co 


o 


+  A/2 
-  A/2 


[JX(P)]  [jY(uvp))  dp, 


for  (-2a.'0  -  A)  <  o'  <  ( -2ajQ  +  A). 


(2.18) 

Combining  F  (co)  and  F„(co)  from  Equations  (2.17)  and 
R  i 

(2.18)  to  give  Fq(co)  =  \  [Fg(o))  +  F,,(co)]  for  quadrature  correlation 

yields : 

-co  +  A/2 

F  (co)  =  Fr(co)  =  ±  /  X(p)  Y(co-P)  dp 

J  -co  -  A/2 
o  7 


co  +  A/2 

1  i  0 

+  ^  X(p)  Y(co-P)  dp, 

J  coq  -  A/2 


for  -A  <  co  <  A, 

=  0,  for  (2ao  -  A)  <  co  <  (2coq  +  A), 

=  0,  for  (-2co  -  A)  <  co  <  (-2co  +  A). 

'  o  '  '  o 


(2.19) 

Thus,  from  Equations  (2.17)  and  (2.19),  F-(co)  and  Fn(co) 

are  identical  for  -A  <  co  <  A  and  F^( co)  is  zero  elsewhere,  while 

F_(co)  contains  second  harmonic  terms  around  co  =  +  2co  as 

K  —  o 
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llluitnt«4  in  Figures  2.1  and  2.5-  It  will  now  be  shewn  that  t  lese 
—coed  harmonic  terns  in  the  regular  correlation  product  spectrum 
produce  sampling  errors  which  are  not  present  in  the  quadrature 
correlation  function. 


B.  Correlation  Accuracy  for  Deterministic  Signals.  In  order 
to  evaluate  the  correlation  function  sampling  errors  when  both  the 
Input  and  references  are  deterministic  signals  (or  t ine- delayed 
versions  of  the  same  deterministic  signal),  consider  first  the  problem 


of  representing  the  integral  of  an  arbitrary  Fourier  transformable 

function  f(t)  by  the  sum  of  uniformly  spaced  samples  f(mT  )  of 

8 

the  function. 

Let  f(t  be  any  Fourier  transformable  deterministic  real¬ 
valued  signal  defined  for  0  <  t  c  I  and  essentially  rero  elsewhere. 
Me  can  represent  f(t)  in  the  time  Interval  T(  by  the  complex 


Fourli 


4  v 


Cn  # 


(2  20) 


where 


c  « 
n 


f(t)  e 


(2.21) 


Next ,  consider  sampling  f(t)  at  a  uniform  sampling  rate  f 


where  I. 


loe  between  samples  in  seconds,  and  averaging  the 


resulting  N  w  samples;  l.e. 
■ 


i 


Since  the  Fourier  series  for  f(t)  in  Equation  (2.20)  is 

_j2«n  (t  i  T)  _j2nn  t 

periodic  with  period  T  and  since  e  =  e 

the  integral  of  Equation  (2.21)  for  the  coefficient  c^  can  he 

rewritten  as: 

T/2  ,2Jtn  . 

~  3  rrt  ^ 

cn  =  J  f(t)  e  dt, 


and  Equation  (2.22)  as 


f(mTs) 


fe  y 


T/2 

j  6 


t 


N  ,2nnm 

Y 


N  ,2nn 

V  JT  “  (1  -  ej2nn 

Substituting  /  e  =  ^ 


0  for  n  /  kN 


ej2nn  r 0  for  n 

7TO - *  ~  j  »  . 

nr  for  n 


(e  ”  -1) 


where  k  is  any  integer,  gives 


I  >  f<"V 


co  T/2  .2«kN  . 

1  \  i  "j  T  ^ 

i  j  f(t)  e  dt.  (2.23) 

k=-oc  -T/2 


Moving  the  k  =  0  term  from  the  right-hand  side  of  Equation  (2.23) 


to  the  left-hand  side  and  noting 


f  gives : 


where  E_  is  the  error  in  representing  the  integral  ^  /  f(t)  dt 


N 


by  the  discrete  sum  ^  ^  f(mTg).  If  the  inverse  Fourier 


m=l 


Transform 

f(t) 

■  J 

F(2«f)  e^2r[tX 

1  df 

is  substituted 

-00 

right-hand  side 

of  Equation  (2.2*0  > 

et 

becomes : 

CO 

r - 1 

rT/2  • 

et 

i 

rn 

y 

/  / 

F(2nf)  df  e  ' 

T 

T 

k=-» 

J  -T/2  j  -oo 

k/0 

00 

V  1 

00 

/ ' 

,T/2  j2nt(f. 

— 

) 

/  F(2itf) 

J- 

T 

e 

k=-oo 

J 

-00 

i 

_  J  -T/2 

k/0 

00 

V 

00 

1 

sinitT(f-kf  ) 

- 

z 

k=-oo 

/  F(2nf) 

-00 

jrT(f-kfs)  df‘ 

j2nkf  t 
s 


dt 


dt 


(2-2: 


k/O 


Equation  (2.25)  can  be  simplified  further  by  writing  F(2itf)  in 
terms  of  its  real  and  imaginary  parts  F(2nf)  =  F(2nf)  +  j.F(2nf) 


2 4 


to  give 


T'Zr.t ) 

r 

k 


sln*T(f-kf  ) 

.?7.kf  r  df 


c 

♦  j  iF(2vf) 


k 


sln*T(f-kf  ) 


The  second  integral  is  zero  because  the  integrand  is  an  odd  function 
of  frequency.  The  first  integral  can  be  reduceo  to  the  Integral  over 
positive  frequency  since  the  Integrand  is  an  even  function  of 
frequency.  Thus,  E^  bee ones 


E 


T 


T(Z*t) 


sln«T(f-kfg) 

“^TTOT7“ 


(*.*«) 


For  a  given  rF(2*f),  T,  and  fg,  Equation  (2.26)  could  be  used  to 
compute  the  error  E^. 

Further  simplification  of  Equation  (2.26)  results  In  aany 
cases  for  signals  with  large  TV  products.  If  ^  «  V  and  rF(2* f) 
is  approximately  constant  for  (fg  -  y)  <  f  <  (fg  »  ^;),  Equation  (2.26) 
approximately  reduces  to 


2  F(2nkf  ) 

r  s 

k/0 


sln«T(f-kf  ) 

3TTST7  df 


If  f 
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the  above  Integral  Is  approximately  equal  to 


1 

T 


for  all  a  C .  Let  t  repreaent  the  iffroilw!  1  _-r.  t  * he  err  r 
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■  a  1  a  1  a  os  os 


•  N  'kT  !»_ ' aT  )  co*  -  kT  sin  u.  aT 
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n  kT  '  n'aT  '  •  •  . (k-a)T  }, 

a  s  nn  a 


r.,  kT  )  K,(aT_ )  *  *  :  ( k-a ) T  ] 

lsls  s 


*1*1  • 


\  *  :(a-k)T  ], 

S2S2  * 


‘1  *  a 


Vz  • 


•\n2 


wtsrs  f  r  k 


\  _ 

r.  s 


*2 


*r>l 


1  2 


Substituting  Equations  (2.38)  and  (2.39)  into  Equation  (2.’*0)  gives 
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The  in-phase  quadrature  correlation  function  for  these  same 
conditions  is,  from  Equation  (2.1*): 


; n ( 2mT  )  y(2mT  )  -  n(2mT  )  $(2mT  )]. 

5  S  5  O 
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Firm  the  above  equations,  4  [(k-m)T  ]  and  4  l ( k-m )T  j  can  be 
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written  as 
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from  Equation  (2.37)  now  becomes: 
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The  mean  value  of  <t>  is  zero  and  the  variance  is 

*1 
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_  M/2  M/2 


2  _  .  2  1 
""  «  =  M2 


-  -2  )  )  [n(2kTg)  y(2kTg)  +^(2kTg)  y(2kTs)] 


k-1  m=l 


[n(2mT  )  y(2mT  )  +  $(2mT  )  $(2mT  )] 
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M/2  M/2 
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s  s  s  s 


+  r?(2kT  )  y(2kT  )  n(2mT  )  y(2mT  )]. 
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(2.42) 


The  first  expectation  in  Equation  (2.42)  is  given  by  the  regular 
correlation  function  result  of  Equations  (2=38)  -  (2.40)  as 


n(2kT  )  y(2kT  )  n(2mT  )  y(2mT  )  = 
s  s  s  s 


4>  M  [2(k-m)T  ]  cos  2cn  (k-m)T 
N1  1  S  °  3 


+  4>N  N  [2(k-m)Tg]  sin  2o)o(k-m)Tg 
1  2 


\  *E  E  [2(k-m)Tg  cos  2coo(k-m)Ts 


+  <>E  E  [2(k-m)Tg]  sin  2^o(k-m)Ts 
1  2 


(2.43)' 
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The  fourth  expectation  n(2kT  )  $(2kT  )  n(2mT  )  y(2mT  )  can  also  be 

s  s  s  s 

shown  to  be  given  by  Equation  (2.1*3).  Similarly,  the  second 
expectation  in  Equation  (2.1*2)  is: 
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s  s  s  s  s  s 
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1  2 


v  [2(k-m)T  ]  sin  2o>  (k-m)T 
s  OS 


+  f>E  E  [2(k-m)Tg]  cos  2u>o(k-m)Ts  >  , 
L  2  I 


(2. kb) 
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and  the  third  expectation  becomes: 


(2kT  )  y(2kT  )  n(2mT  )  y(2mT  )  =  [n'(2kT  )  n(2mT  )] 
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cam  now  be  obtained  from  Equations  (2.1*2)  -  (2.1*5 )» 
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2 

A  comparison  of  oR  from  Equation  (2.1*1)  for  the  regular 

2 

correlation  function  and  from  Equation  (2.1*6)  for  the  quadrature 

correlation  function  shows  the  presence  of  second  harmonic  terms  in 

2  2  2 
aR  which  do  not  appear  in  .  Thus,  while  depends  only 

upon  the  relation  between  the  sampling  period  T  and  the  signal 

s 

2 

bandwidth  W,  aD  also  depends  upon  the  center  frequency  u)  of 

the  signal  spectrum.  This  is  illustrated  in  the  example  derived  in 

the  Appendix  for  the  case  where  n(t)  is  bandlimited  white  noise 

with  a  flat  power  spectrum  for  (fQ  -  -|)  <  f  <  (fQ  -1-  |[) ;  For  this 
2  2 

example,  aD  from  Equation  (2.1*1)  and  cr»  from  Equation  (2.1*6) 

K  ^  f 

are  plotted  in  Figures  2.6  and  2.7,  respectively,  for 


-  9  and 


a  21V  product  of  Figure  2.6  shows  that  the  regular  correlation 

variance  can  be  as  much  as  2  to  3  IB  higner  than  the  quadrature 
correlation  variance  of  Figure  2  7  depending  -upon  the  relationship  of 


the  sampling  rate  to  fm  and  W.  For  this  example,  not  only  is  the 
effect  of  the  second  harmonic  missing  in  Figure  2  ",  but  the  quadrature 
correlation  function  variance  is  seen  to  be  constant  at  about 


2 


for  all  f  >  W. 

s  — 

To  summarize,  for  the  quadrature  correlation  function,  the 

sampling  error  E  was  found  to  be  zero  and  the  correlation  variance 

2  a  A 

oQ  was  found  to  be  a  constant  independent  of  f  for  f  >  A/2n. 

However,  for  the  regular  correlation  function,  both  the  sampling 

error  and  the  correlation  variance  were  found  to  have  undesirable 

second  harmonic  terms  at  some  values  of  f  . 

s 


2.3  Comparison  of  Quadrature  Correlation  Functions  for  Clipped 
and  Unclipped  Signals 

In  this  section,  we  wish  to  consider  the  properties  of  the 
quadrature  correlation  function  for  a  special  class  of  functions- 
clipped  narrowband  signals.  The  relationship  between  a  real-valued 
function  x(t)  before  clipping  and  its  value  u(t)  after  clipping 
is  given  by 


u ( t )  =  sgn  x(t) 


I  1  for  x(t)  >  0 
]  -1  for  x(t)  <  0. 


i*2 


Hard  limiting  thus  destroys  all  amplitude"  information  in  x't) 
and  retains  the  phase  "  information  in  the  sense  that  only  the  tero 
axis  crossings  of  x(t  are  preserved  in  u(t). 

The  analytic  model  for  the  quadrature  correlation  function  for 
clipped  signals  shown  in  Figure  2.8  differs  from  the  quadrature 
correlation  function  model  for  unclipped  signals  of  Figure  2.2  only 
by  the  fact  that  u(t)  and  v(t)  are  hard  limited  versions  of  the 
signals  x't)  and  y(t).  The  quadrature  correlation  function, 
Equations  '2  1 )  -  [2  5)>  can  thus  be  used  to  describe  the  clipped 
correlation  functions  of  Figure  2.6  by  replacing  x(t)  by  sgn  x(t) 
and  y't )  by  sgn  y't) . 

We  wish  to  compare  the  clipped  and  unclipped  quadrature 

correlation  functions  of  x(t)  and  y(t)  for  two  cases  : 

]1]  x't)  =  o  s(t)  and  y(t)  -  s(t  -  t)  are  both  equal  to  the 
s 

same  narrowband  gaussian  signal  s(t)  but  with  different  relative 

amplitude  a  ,  and  2]  y(t)  -  s(t)  and  x(t)  -  c  n(t),  a  member 
s  n 

of  a  wide-sense  stationary  gaussian  random  process.  Let 

n(t)  N,  (t)  cos  u.  t  -  IL't)  sin  o  t 
1  o  2  o 

and 

s(t)  -  E. (t)  cos  u  t  -  E_ (t )  sin  u  t, 
v  1  o  2  o 

where 

NL2(t)  -  N22(t)  -  1 
EL2(t)  -  E22(t)  =  1 


<  ► 


and 


E  ft)  K  (t  -  t)  =  0  for  i ,  j  =  1,  2. 

The  values  c  and  c  account  for  differences  in  relative  amplitude 
s  n 

between  x(t)  and  y(t). 

Let  the  quadrature  sampling  rate  be  at  the  Shannon  rate  of 
f  =  -—r-  =  W  sample  pairs,  x(2nT  ),  x(2nT  )],  per  second. 

S  Cm  X  S  S 

s 

The  Hilbert  transform  of  s(t)  is  given  by 

s(t)  E, (t)  sin  u  t  +  E„(t)  cos  u  t 
1  o  c  o 

and  the  pre-envelope  of  s(t)  by 
a 

s(t)  -r  js(t)  -  E.,(t)  1  jE2(t)l.cos  uQt  +  j  sin  u^t]  =  f(t)  e 

where  f(t)  is  the  low  pass  modulation  function  whose  spectrum  is 
equal  to  the  positive  spectrum  of  s(t)  translated  down  to  zero 
frequency  and 

*ff(0  -  f*(t)  f(t  -  t)  -  2:E1(t)  Ex(t  -  t)]<-  J2fE1(t)  E2(t  -  x)] 


2  s 


E1E1 


(t) 


2  ®  p  y  (  0  . 

E1E2 


Consider  first  the  quadrature  correlation  function  for  the 


unclipped  case  defined  by 


<o 


M/2 

*q(Mts )  -  i  X  ;x(2aTg)  y(2nTs)  +  x(2bTs)  y(2nTs)]f 

o=l 

M/2 

(MT)  =  i  1x(2bT  )  y(2mT  )  -  x(2mT  )  y(2nT  )], 

S  n  5  5  S  5 

m-1 

and 

w  =  v1".1  *  -  V".’- 

Considering  Case  1,  if  we  let  x(t)  =  o  s(t)  and  y(t)  =  s(t  -  t 

s 

Zq(t)  will  be  the  auto-correlation  function  of  s(t).  Computing 
z^(x)  at  increments  of  t  equal  to  the  spacing  between  samples, 

2Ts,  gives 

M/2 

zQ(2nTs)  "  E  x(2mTs)  -  jx(2BTs)][y(2mTs)  +  jy(2mTg)J 

m=l 

M/2 

--  ^  )  ts(2mTs)  -  js(2mTs)]'s(2oTs  -  2nTj 

m=l 

♦  js(2mT  -  2nT  ) ] 

5  o 

J2w  nT  'P 

■iv  m  \  f*(2nv f(2BTs •  2nV 

m=l 

a  j2oj  nT 

=  re  °  5  zff(2nTs}»  (2>ll8) 


where,  by  Equations  (2.27)  and  (2.28),  z^(2nTs)  is  equal  to 
1  T 

i  f*(tx  •  f(t  -  i)dt  for  f  >  W  and  t  =  2nT  .  Taking  the 

1  S  “  s 

o 

real  and  imaginary  parts  of  Equation  ( 2 .1*8 ) : 


:q(2"V  -■  -y'e,f.(2"ts)  c°s  *■„*'  •  v,(2nv  5in  ’VV- 


and 


(2.«»9) 


V2nV  *s  ::E1E1f2nTs)  Sin  2;‘onTs  +  JE1E2(2nTs)  COS  ^onTa]- 


(2.50) 


The  correlation  function  envelope  is  given  by: 


*V*>V  Y=Q2(2nTs)  * 


=a(2nT,)| 


V  (2nTJ* 

fcl  2 


(2.51) 


The  maximum  value  of  b^(2nTs)  occurs  at  n  -  0  and  is  given  by: 


^Q^max  “  W°;  *  V 


(2.52) 


For  the  clipped  case,  the  auto-correlation  function  of  s(t) 
can  be  found  by  using  the  relationship  37]  between  the  cross¬ 
correlation  function  of  two  Gaussian  signals  before  and  after  clipping: 


VvN\ 

e.b(,) 


I  sln'1  °«b  ( '  > 


(2.53) 


where 


croeo-correlat Ion  function  between  a(t)  and  b(t)  after 
clipping, 


1*8 


» 


a* 

i  ► 


(£,L  '  \(0)  •  *  1  (2.56) 


Equation!  (2.56)  and  (2.52)  lllut trait  a  baaie  noraalltatlon 
proparty  of  clipping;  l.a.,  tha  correlation  function  for  clipped 
signal*  ha*  a  rang*  of  -1  to  +1,  while  for  the  undipped  signal,  the 
range  io  -<j  to  ,  a  function  of  the  signal  amplitude. 

For  Case  (2),  if  y(t)  *  s(t)  and  x(t)  *>  onn(t),  a  member 
of  a  stationary  gausslan  random  process,  the  average  value  of  the 
quadrature  correlation  function  for  the  unclipped  case  Is  tero  and 
Its  variance  as  given  by  Equation  (2.1*6)  Is: 


»„2 

TT 


2 

M 


m*  1 


(• 


NjN^(^aTs ) 


F  ( 2mT  ) 
T1  8 


(2oT8)  #e  e  (2oTg))(l  - 


2m. 

rK 


(2.57) 


If  n(t)  Is  assumed  to  have  a  rectangular  spectrum  of  width  W  Hz, 

/\ 

Figure  2.7  shows  that,  for  any  f  >  W, 

0  * 


(2.58) 


where  T  is  the  signal  duration  in  seconds. 

The  correlation  function  variance  for  the  clipped  signal  case 

MV. 

2 

can  be  found  by  computing  cQ  from  Equation  (2.1*2)  as  a  function  of 

O 

u(t)  and  v(t)  and  then  using  Equation  (2.53)  to  express  0^  in 
terms  of  s(t)  and  n(t). 
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M/2 


M/2 

y 

m*l 


/W\  AV\  . 

♦  ♦u^(2(k-m)Ti]  *vA(2{k.m)TB)} 


2m  ^ 

(l  '  r)l*uu(2mTn)  \v(2mV 

m  1 


/wv  /wv 

4  4u(J(2mT«)  >vO(2nTB)| 


2m » 

r> 


» 

nn 


( 2mT  ) 
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Jlsin"1 


$ 

BE 


(2mT  )) 

O 


*  *(2nT  )  .  >* 

♦  .'ain"1  nn  2  8  ][ain  %g(2mTa)H.  (2.59) 

’n  * 

A  plot  of  the  variance  given  by  Equation  (2  59)  for  a  random  process 

n(t)  with  a  rectangular  spectrum  of  bandwidth  W  Hz  and  center 

f 

frequency  tQ  Hz  is  shown  in  Figure  2.9  for  ^  9  and 

M  -  21V  =  50  The  equations  for  this  example  are  derived  in  the 

Appendix.  Figure  2.9  shows  that,  for  the  Shannon  sampling  rate  of 
A 

f  =  W,  the  variance  is  given  by 

8 

MW 


2  _1_ 

°Q  "  21V  ’ 


(2.60) 


but  if  the  sampling  rate  is  increased  sufficiently,  a  variance  of 


approximately 


MM 

2 

°Q 


1 

1*TV 


(2.61) 


Clipped  Signal  Quadrature  Correlation  Function 
Variance  Vs.  Sampling  Rate 


can  be  achieved.  The  small  ripples  in  Figure  2.9  occur  when  integer 
multiples  of  the  sampling  rate  are  equal  to  the  second  harmonic  of 
the  signal  center  frequency.  Because  of  the  non-linearity  of  the 
limiters,  these  ripples  were  not  completely  eliminated  as  were  those 
of  Figure  2.7  for  the  unclipped  signal  case. 

For  comparison  with  Figure  2.9,  the  variance  for  the  regular 

clipped  correlation  function  is  shown  in  Figure  2.10  for  the  same 

example.  The  curve  of  Figure  2.10  shows  the  same  asymptote  for 

large  f  as  the  quadrature  correlation  curve  of  Figure  2.9>  but  it 
*  8 

exhibits  ripples  of  about  3  dB  instead  of  approximately  0.5  dB 
for  the  latter  case. 

A  comparison  of  Figure  2.7  and  2.9  shows  that  clipping 
permits  a  reduced  correlation  variance  for  random  signals  if  sampling 
is  increased  sufficiently  beyond  the  Shannon  rate.  In  summary,  the 
two  principal  effects  of  clipping  are  seen  to  be  correlation  function 
normalization  and  reduced  correlation  variance  for  high  sampling 
rates . 
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CHAPTER  III 

QUADRATURE  CORRELATION  FUNCTION  PROPERTIES  FOR 
SINUSOIDAL  SIGNALS  AND  UNIFORM  SAMPLING 

3.1  Quadrature  Correlation  Function  for  Sinusoidal  Signals 

In  this  chapter,  the  uniform  sampling  quadrature  correlation 
results  of  Chapter  II  for  narrowband  functions  with  large  TW  products 
will  be  extended  to  the  case  of  sinusoidal  signals.  After  examining 
the  correlation  properties  of  sinusoidal  signals,  the  effects  of 
uniform  sampling  and  hard  limiting  will  be  derived  separately  and, 
finally,  the  combined  effects  of  both  sampling  and  hard  limiting  will 
be  examined.  The  problem  of  spurious  odd  harmonic  correlation 
resulting  from  clipping  and  sampling  will  be  investigated  in  detail 
and  a  solution  to  the  problem  by  the  use  of  non-uniform  sampling  will 
be  introduced. 

Consider  the  quadrature  correlation  function  model  with  finite 
integration  time  T  shown  in  Figure  3.1  for  two  infinitely  long 
sinusoidal  signals  with  different  frequencies.  Let  the  input  and 
reference  signals  be  given  by 


and 


x(t)  =  COS  CD^t,  -oo  <  t  <  oo  , 


y(t )  =  cos  co  t, 


-00  <  t  <  00, 
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where  u>,  -  2«f,  and  ta  -  2nf  are  the:  sinusoidal  frequencies. 

110  0 

As  discussed  in  Chapter  II,  the  real  and  imaginary  parts  of  the 
complex  correlation  function  z^(t)  of  Figure  3.1  represent  the 
"in-phase"  and  "quadrature  phase"  correlation  components  and  z^t) 
represents  the  envelope  of  the  correlation  function.  From  Figure  3.1, 
z1(t)  is: 


2i(t)  *  h  .. 


T/2  jcujt  -ju>  (t-i) 
e  e 

-T/2 


dt 


L  Jo)qt  sin  it  T^-f^) 

2  6  x  T(frfo} 


The  correlation  envelope  peak  at  t  -  0, 


ai(°> 


sin  n  T(f1-fQ) 

rt  T( f  -f  ) 
v  1  o 


(3.2) 


is  plotted  in  Figure  3.2  as  a  function  of  the  frequency  difference 

(f^-fQ).  This  plot  shows  that,  except  for  frequencies  f^  within 

about  -jj  of  the  fixed  frequency  fQ,  the  correlation-envelope 

peaks,  as  a  function  of  input  frequency  f^,  never  exceed  the 

s  i  n 

sidelobe  peaks  of  the  -  function.  Only  the  first  two  peaks 

X 

are  greater  than  10  percent  of  the  main  lobe  value  at  f^  =  fQ. 

For  large  T,  therefore,  this  correlation  function  has  a  very 
selective  frequency  characteristic,  since  it  has  appreciable  amplitude 
only  for  frequencies  f^  which  are  within  about  -jj  of  the  reference 
frequency  f  . 
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The  correlation  function  given  by  Equation  (3.1)  assumes  an 
input  signal,  cos  o^t,  of  infinite  duration.  For  a  sinusoid  signal 
of  finite  duration  T,  an  additional  linear  weighting  factor  is 
present  giving 

1  .J'V  sin  ”  T<fl-fo>  „  111  , 

Zjl'O  2  e  n  T(f.-f  )  "  T  '  • 

1  o 

Note  thu*-  the  correlation  peak  at  x  0  is  the  same  as  that  given 
by  EquatiCs  (3.2)  for  the  infinite  duration  case,  justifying  the  use 
cf  infinite  duration  sinusoids  for  finding  the  frequency  response 
for  small  t  sucn  at  t«T. 

3.2  Uniform  Sampling  Effects 

Next,  consider  the  quadrature  correlation  of  sinusoidal  signals 
with  uniform  sampling.  Figure  3.3  shows  a  quadrature  c  rrelall  >n 
function  model  with  uniform  sampu  g  intervals  f  I’l  sr<  <4.' 
with  input  and  reference  sinusoidal  fi.-mals  f  «  c  *  ,  r*-  p*-»iv»i, 
The  complex  quadrature  correlation  funrt .  ■  '  f  »  infnu'r 

duration  Input  sinusoids  Is: 
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*2(t)  = 


1 

M 


j^l2mTs  -ju>  asf  0%T 

e  e  e 


0=1 


*V 


m-1 


eJV  tin  n  MT  (f.-f  )  jT  (f  +  1 )  2«(f  -f  ) 

ir  Tirrmr^n  *  (3-3) 

•  l  o 


,M 


where  -  Is  the  number  of  staple  ptirs  taken  during  the  finite 
averaging  time  of  T  seconds.  The  correlation  function  for  i  -  0, 


*2(0) 


1 

I 


sin  «  T(f,.f  ) 
i  o 

£  M  sin  «2T#(frfo) 


(3.1*) 


shown  in  Figure  3 .*>  approximates  the  *■■  ■*  envelope  characteristic 


of  Figure  3.2  but  with  a  repetition  Interval  of 


? 


1 

2 r  •  M 

s 


expected  from  sampling  theory  36].  The  sldelobe  peaks  midway 
between  the  melnlobes  are  equal  to  about  1  /X,  where  M  is  the 
total  number  of  samples  This  correlation  function  also  has  a  very 
selective  frequency  characteristic  as  long  as  the  Input  frequency 
f^  does  not  differ  from  the  reference  frequency  fQ  by  greater 


than 


f 


Herts . 
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3.3  Generation  of  Odd  Harmonics  by  Clipping 

Next,  consider  the  effect  of  hard- limiting  on  the  correlation 
of  two  sinusoidal  functions  as  shown  in  Figure  3.5.  The  clipped 
sinusoidal  signals  x(t)  and  y(t)  are  square  waves  of  unit 
amplitude  which  can  be  represented  by  a  Fourier  series  of  odd 
harmonics  of  the  fundamental  input  frequencies  f  and  f^, 


x(t)  =  sgn(cos  o^t) 


and 


y(t)  =■  sgn(cos  u^t) 


h 

n 


k 

n 


l 


k=0 


{iwT)  c°s  (2n  +  1>  V' 

T|ji^  cos  (2k+l)  V. 


If  the  boxes  labeled  "Hf '  of  Figure  3.5  are  assumed  to  shift  the 
input  square  wave  by  90  degrees  of  the  fundamental  frequencies 


“L 


and  to  ,  then  the  pre- envelopes  p  (t)  and  p  *(t-x)  become: 


px(t) 


x(t)  + 


A  ji  V  [  -nn  j(-l)n(2n+l)  cu  t 


(3.5) 


n-0 


and 

Py*(t-T) 


y(t-T)  -  jy(t-T) 

i  V  e'J 

it  /,  (2kH) 

k-0 

(-l)k(2k+l)  coo(t-r) 

> 

(3.6) 
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and  the  complex  correlation  function  is: 


Zj(t)  -  2T 

T/2 

il/2 

px(t)  p  *(t-T)  dt 

It 

00 

Y 

n=0 

«  j(-l)k(2k*l)  . 

/  ('l)  '  |2n>l)(2kt!) 

k=0 

(t 

T/2 

-T/2 

e  dtjt 

where 

(3.7) 


-  [  (-1 )°( 2n*l )  w ^  -(-l)k  (2k*l)  wq)  represents  differences 

between  the  odd  harmonics  of  and  w  for  given  values  of  n 

and  k.  Table  3.1,  which  shows  som  typical  values  that  can 

as  sum,  contains  all  the  values  of  P.  with  amplitude  weighting 

an 

fMtor**  rsCTTTarn  •  «r~i#r  xtmti  °-1- 

Mote  from  the  definition  of  that  the  odd  harmonic 

frequencies  alternate  sign;  i.e.,  Is  positive,  is  negative. 

Is  positive,  etc.  This  Is  the  result  of  the  fact  that  the 
boxes  In  Figure  3.?  tlee  shift  s(t)  and  fit.)  by  a  tine  equivalent 
to  90  degrees  of  their  fundaeental  frequencies  •  and  .  n>ve»er. 
this  sane  t  1m  shift  Is  equivalent  to  2  TO  degrees  (or  .90  degrees) 
for  the  third  hereonlc  exponent  of  the  square  wave,  to  135  degrees 
or  *90  degrees)  for  the  fifth  harmonic  empnent,  etc.  Thus, 
positive  Hilbert  transferee  are  obtained  for  the  fundaeer.tal ,  fifth 
hereonlc ,  ninth  harmonic ,  etc  ,  while  negative  Hilbert  tranaf  res 


id 


(■ 
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TABLE  3.1 

ODD  HARMONIC  DIFFERENCE  FREQIXNCIES  WITH 
WEIGHTING  FACTORS  OYER  10  PERCENT 
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are  obtained  for  the  third  harmonic,  seventh  harmonic,  etc.  This 
phenomenon  appears  as  alternating  positive  and  negative  odd  harmonic 
frequencies  in  the  spectruns  shown  in  figure  3.6  of  the  pre-envelopes 
Px(t)  and  Py*(t)  obtained  by  taking  the  Fourier  transforms  of 
Equations  (3.5)  and  (3.6)  end  given  by: 

m 

V">  o  p*(l)1  ‘  t  '('l)n  (2n#l)  W1] 

n^O 


and 


p  •(-)  *  £/ 


(-1  )*  (2k*  1  )  -  ). 

0 


her  forming  the  Integration  indicated  ir.  Equation  3.7)  (ires 
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00  oo 


13(t) 


.  s,  v  y  <■» 


j(-l)k(2k+l)  w  t  Sin  2 


n+k  e 


( 2n+l )  (2k+l) 


n=0  k=0 


sin(w^  -  “0)  g  jV 

L  (-x-“o>I  ' 


—  .  „  T 

^kn  2 

input  reference 
harmonic  harmonic 


where  FUND. 


FUND. 


sin^  +  «,)  2 

f - * 

(**>1  +  “o*  2 


third  fund. 


!  *in(wl  *  5wo)  2  _‘J3V 

-  T  T  e 

(wl  *  3uo)  2 


FUND.  THIRD 


T  -  wo)  2  JV 

♦  T  - -  • 

’  (5^  -  «"0>  2 


FIFTH  FUND. 


ln(u,  -  >•»  )  T 


l  >Wwl  *  ^o;  2 


FUND  FIFTH 


(3.8) 

where  the  individual  c exit rl but l one  of  all  combination*  of  Input  and 
reference  hamonic.  can  be  eeer.  For  the  .pedal  cue  of  ^  -  wq 

and  ,  -  0.  all  the  tern,  in  Equation  (1.6)  are  nefllfible  eacept 

when  k  *  n.  tlvlnf: 
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oo 

m°>  =  %  y 


"  „T0 


-  =  1  for  u)q  =  uv 


Likewise,  for  -  -3^  and  t  =  0,  all  terms  are  negligible 

except  when  n  -  (3k+l),  giving: 


.(0)  -  h  t  s 


n2  3  (2n+l ) 


-1  -1  -  , 

— 2  -  -  f°r  “o  *  -5V 


If  we  consider  just  the  peaks  of  zg(0)  ior  frequencies  at  which 
=  0,  we  can  express  the  contribution  due  to  the  principal 

pairs  of  input  and  reference  odd  harmonics  by: 


*2«» 


peak 


V  TsSnT  *  K  •  (-»>-  «nn)  V 

n-0 


}  jferfj  n«*»0  -  (-*>k  (2k+])  (3-9) 

k  1 


as  shown  in  Figure  3.7.  Equation  (3.9)  neglec.s  differences  between 
higher  order  odd  harmonics  of  u»o  and  tie  most  significant  of 

which  (for  n  -  3,  k  -  5)  has  an  amplitude  of  .  niy  ^  . 

As  shown  in  Equation  (3.9),  the  principal  terms  are  of  two  types: 

(1)  differences  between  the  input  fundamental  w.  and  the  odd 
harmonics  of  (the  terns  indexed  by  n),  and  [?)  differences 
between  the  reference  fundamental  u  and  odd  harmonics  of  u  (the 
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Figure  3.7  Frequency  Selectivity  of  Correction  Functlc n  for 
Clipped  Sinuaoldal  Slgnela 
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1.#.,  to  1/!  for  third  terwfiici,  1  /*>  for  flflfe  temxte*,  etc 


The  nest  section  will  consider  the  effect  of  uniform  uall 


on  ihe  modal  shown  in  flfur*  3.3.  It  w.  11  be  itojwn  that  for  the 


1*4  (ui (  t he  o44  totmmle  larva  shown  in  Plfurs  5.7  *t 


-1)  !2n'l  -  will  be  repealed  at  Multiples  of  the  seep  l  Inf 


frequency  2»f( ,  end  the  sub -harmonic  teres  vt 
will  be  repeated  vt  Multiples  of 


5.1  Effect  of  tolform 


L.'_  1H?  ■fill  1M1L7LHI 1 


To  esse  Ins  the  effect  of  uniform  snMpllr*  of  clipped  sinusoids , 


consider  the  quadrature  correletlon  function  Model  shown  la  Flpure  3.6 


The  express  1  jns  for  p  ( t  )  and  p  *(t.  *  '  ere  clean  by  Equations  (3. 9  ) 

®  y 

and  (3.6).  The  sanpler  output  signals  are 


p  f  2m  T  ) 

*S  S 


j  -1  2n*  1 


(5  10) 


p  •(  2m  T  -  t  ) 

y  • 


•“  ,  ,  >h  • ; ! • 1 )  (2h*l )  -  (2m?  .* ) 

T  •  ' 


(5.11  ) 


where  2T 


l /fr  Is  the  quadrature  pairs  asepllrq  Interval.  The 


complex  correlation  function  is 
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H(T)  =  i  )  Px<a»Ts)  P  *(2»TS  .  T) 


00  00 


tl(-l)K(2k+l)  u>  t 


16  \ 

T  / 


)  (-1)**  Z- 


(2n+l)  (Lk+i: 


n-0  k-0 


T  M^2  i2»T  0 

1  \  j2“Vkn 

•  H  /  e 

m- 1 


(3.12) 


where 


(-1  )n(2n*  1 )  wx  -<-l)k(2kH)  wqJ  4*  before. 


Putting  the  bracketed  series  of  Equation*  (3.12)  in  closed 


for*  giver 


w 


,n*k  e 


JT.  <5  *  11  V 


n  0  k  0 


’  sin  |  W|  flta  1  j(.nk(2k*i)  ->ot 

■JT7^T7“J 


(3.13) 


where  HT 


' 
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r 


Writing  Equation  (3.13)  in  terms  of  individual  odd  harmonics 


contributions  gives: 


,  ,  8  r  sin(ui  -  uo>  t/s 

2  |_M/2  sin(u1  -  u  ) 


JT  (|  +  l)(u  -  u  )  I  T 

e  e 


l  sin(3ui1  +  uq)  T/2  -jTg(|  +  l)(3w1  +  u)q) 

3  M/2  sin(3u1  +  u  )  Tg  e  * 


1  sin(u1  +  3ujq)  T/2  jTg(|  +  L^U)i  +  3u,0)  -j5w0T 

3  M/2  sin(w1  ♦  3uq)  Tg  *  * 


L  sin(5w1  -  l>o)  T/2  jTg(|  +  lM^  -  wQ)  jujQi 

+  5  M/2  sin  (5^  -  <*>Q)  Tg  e  e 


L  sin(u)1  -  5u>q)  T/2  jTg(|  *  1)^  -  VD)  J%>0i 

+  5  M/2  sin(u  -  %  )  t  *  * 

1  0  8 


(3.14) 

This  equation  shows  that  *^(0  contains  repetitive  terns  due  to 
sampling  similar  to  those  shown  in  Figure  3.4  about  each  odd  harmonic 
difference  frequency  satisfying  the  relation  0. 

Neglecting  terms  with  amplitudes  of  the  order  of  1/1}  and 
sattller,  an  equation  analogous  to  Equation  (3  9)  can  be  written 
approximating  principal  peaks  of  (0 )  as 


tk 


\(o)) 


peak 


I 


2k*l ) 


>  [uj  -  (-ir(2k«l)  wc 


a 

“V 


k-0 


n-  1 


P  u 

0 


( -1  )n(  2n*l ) 


(S  15) 


If  the  reference  frequency  w  -  2*f  In  Equation  (3.15)  1*  held 
fixed  while  the  input  frequency  2«f^  la  varied,  the  complex 

pattern  of  odd  harmonic  correlation  tema  typified  by  thoee  abown  in 
figure  3.9  for  the  f undaaental ,  third  and  fifth  harmonica  will  result . 
Thia  figure  llluatratea  that,  becauae  of  uniform  sampling,  aany  odd 
harmonic  correlation  peaks  are  produced  for  input  frequencies  f 
which  differ  from  f  by  only  a  fraction  of  the  aampllr^  rate  . 

This  could  be  a  highly  undesirable  situation  If  a  correlation  function 
with  good  frequency  selectivity  la  necessary. 

The  exact  location  of  the  Input  frequenr les  f  that  give 
rise  to  spurious  odd  harmoc  lc  correlation  peaks  as  a  function  of 

A 

f  and  f(  la  shown  In  figure  3.10  Thia  figure  shea  that  In  any 

of  the  fixed  fl  frequency  Intervals  (h  •  1)  r  of  lee^th 

where  h  is  an  arbitrary  integer,  t  re  re  are  iw  fundamental, 

4  third  harmonic,  6  fifth  harmonic  ,  etc.,  frequencies  giving  spurious 

co-relatl  ns  for  flaed  values  of  f  and  ?  ,  the  l  vati  a  of  eech 

os 

of  these  response  frequencies  can  be  located  from  figure  3. 1C  once  the 


petition  of  the  fintettlAl  rviparjt  f  ♦  Kf  '  for  tea*  Integer 

M  it  eatafellthed  and  It j  dlttanet  ■  fra  the  left-hand  add*  -f  the 
Interval  if  It  tettrvlnH.  TSa  aifnlflcance  of  ttt  al taraat  1  n* 
positive  and  Mfitln  odd  harmonic  fr»qjtnel«i  ih»r.  in  Flfure  3  6 
It  apparent  la  FI4  xa  3. 1C  la  that,  If  the  f  rifleman'  a  1  rtipatt 
frtfutnt/  «tn  t  rt4  to  the  npht  tjr  ta  auat  It,  tM  third, 
seventh,  tlmntfc,  etc . ,  harmonic  frt^utwtti  would  art  to  the  left, 
•hilt  tht  fifth,  ninth,  ttc.,  harmonic  frepuenclea  «ru  1  d  move  to  the 
rljpt  fc/  tht  taunt  indicated  ia  Fifure  3  10 


3 .5  C.rrtlHUn  function  <ar  lance  for  jarr’^dtnl  Kvdt  Prxtn 

Tht  correlation  function  variance  0  ‘ ,  defined  at  the 

% 

rnlawt  of  tit  real  part  f  * )  of  Figure  3.S  *ttn  tht  input 

■  (t  n  t  la  a  f  1  Ipped  narr  wtanl  rand  a  pr«eaa  and  /  t  '  la 

a  clipped  ainutoldal  function  la  <lvtn  fcj  Lputti  n  aa 
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bandwidth,  the  variance  of  figure  3.11  approaches  the  f  st  sampling 
rate  asy^totic  value  of  l/kTV  sore  slowly  than  the  corresponding 
large  IV  case  of  figure  2.9- 

-  -6  Use  of  Xonuniform  Saa pltng  to  Eliminate  Odd  Harmonic 

Wespcnses 

In  the  nest  two  chapters,  we  will  consider  the  use  o" 
nonunlfora  t  1m  s sapling  to  reduce  the  effect  of  spurious  o  i  harmonic 
correlations  caused  by  the  clipping  and  sampling  of  sinusoid.  .  The 
general  sampling  technique  eons.dered  is  that  of  “Block  Sampli-V'* 
l.e.,  the  use  of  short  sequences  (or  blocks)  of  nonunifone  tiao 
»a*ple*  which  are  repeated  at  periodic  Intervals.  The  nonuni fjrm 
stapling  destroys  the  synchronise  between  the  sampling  periodic.' ty 
and  the  signal  pe-.odiclty  which  cause  the  odd  harmonic  responses  for 
the  uniform  sampling  case.  Some  general  theoretical  properties  >.f 
block  sampling  will  be  developed  in  Chapter  IV,  and  the  use  of  b  lock 
sampling  to  reduce  spurious  odd  harmonic  effects  in  quadratuie 
rorrelatlon  functions  of  clipped  sinusoidal  signals  will  be  discur Jed 
in  Chapter  V 


CHAPTER  IV 


SPECTRAL  ANALYSIS  )E  NONUNIFORM  BLOCK  SAMPLING 

L . 1  I ntroduct I  ■  i . 

The  line  nf  repetitive  nonuniform  sampling  (blxrk  s  amp  1  i  ng 1 
has  been  considered  for  various  applications  in  the  literature. 

Tou  [36]  uses  Z-transforms  to  analyze  sampled-data  control  systems 
with  "cyclic  variable-rate"  sampling.  Yen  rLl]  and  Kohlenberg  ' 
derive  sampling  theorems  for  the  reconstruction  of  signals  which 
have  been  sampled  by  various  nonuniform  sampling  schemes,  one  of 
which  is  "recurrent  nonuniform  sampling."  In  this  chapter,  we  will 
examine  some  of  the  basic  properties  of  block  sampling,  concentrating 
on  the  sampling  frequency  spectrum  properties.  The  importance  of 
studying  the  sampling  function  spectrum  is  evident  from  the  fact  that 
the  spectrum  of  a  sampled  signal  is  equal  to  the  convolution  of  the 
spectrum  of  the  unsampled  signal  with  the  sampling  function  spectrum. 

In  Chapter  V,  the  sampling  spectrum  results  from  this  chapter  will  be 
applied  to  the  problem  of  reducing  odd  harmonic  correlation  peaks  fur 
a  correlation  function  between  two  sinusoidal  signals. 

Before  deriving  the  block  sampling  time  and  frequency  equations, 
the  analogous  uniform  sampling  equations  will  be  obtained  for 
comparison.  Consider  the  infinite  train  of  uniformly  spaced  ideal 
sample  pulses  shown  in  Figure  h.l  and  given  by  the  sampling  function 
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where  T  /i  la  the  average  Interval  taden  »tr  ne  t  1  ca  f 
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length  T  .  Since  (  t  '  la  a  perl  41c  functl  m  if  perl  4  7^,  It 

can  oe  represented  If  lta  Purler  series  e spans  1  m  t/ 
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A  corollary  to  Equal  1  on  (A. 14'  la: 


Si. a 


for  all  k. 


Mtii'fa  follows  directly  tram  c~mp\.*x  conjugate  pro  part  its .  Tht 
• /Mr try  pr^crtln  of  Equation*  4.14)  and  (4.9)  art  Illustrated 
1b  Figure  4 .} . 

An  araal nation  of  tha  periodicity  condition  of  Equation  (4.12) 

will  llluatrata  tha  ralatlonahlpa  between  tha  nonunifor*  sampling 

times  ^t  )  *  and  tha  aaapling  spectru*  coafflclanta ,  a^. 
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THEOREM:  Tha  function  a.  -  jj  a  P  is  periodic  If, 
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PROOF: 
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P 

n  and  represent  these  rational  numbers  as  the  ratio  of 
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two  integers,  i.a.,  let  y  =  ,  a-  - 

P  1  P  2 

tN  Sj, 

y-  y  where  a^  and  b^  are  integers  for 

P  N 

1  =  1,  2,  ...  ,  N.  We  need  only  exhibit  one  value  of  K* 


•  i 


(not  necessarily  tha  alnimua  K  defined  by  Equation  (4.10)] 
which  will  satisfy  Equation  (4.12)  for  all  n  to  prove 
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periodicity.  One  such  K*  is  K'  =  •  bg  . . . 

N 

=  n  b. .  To  show  necessity,  consider  Equation  (U.ll). 
i=l  1 
t 

If  is  an  irrational  number  for  any  n,  then 
P 

t 

-J2«K^r 

e  ^  could  not  be  unity  for  that  n  and  periodicity 
of  would  not  be  possible. 


In  practice,  where  sampling  times  cannot  be  determined  with  infinite 
accuracy,  approximate  periodicity  rather  than  precise  periodicity 
must  be  considered.  This  topic,  which  could  be  an  area  of  further 
study,  is  not  considered  in  this  thesis. 

Assume  that  a  set  of  sampling  times,  £tn  ,  is 

specified  and  that  we  wish  to  determine  its  spectral  period  K.  Let 
K'  be  the  smallest  common  demoninator  simultaneously  satisfying  the 
N  equations : 


T  K' 
P 


I 


for  some  set  of  integers  Nh ,  i  -  1,  2,  N.  Then,  the  spectral 

period  K  is  equal  to  K' .  This  is  equivalent  to  saying  that  K  is 

the  smallest  number  of  equally  spaced  time  instants  into  which  each 

block  sampling  interval  T^  can  be  divided  such  that  each  of  the  N 

T  "\  N 

nonuniform  sample  points  ^  t  j  ^  occurs  at  one  of  these  K 
uniformly  spaced  time  instants.  This  is  illustrated  in  Figure  4.6. 

If,  instead  of  the  sample  times  On)  being  specified,  a 
desired  periodicity  K  is  given,  the  acceptable  sets  of  sampling 
times  giving  this  periodicity  can  be  found  by  dividing  the  block 
sampling  interval  into  K  equally  spaced  time  instants  T^/K 

seconds  apart  as  shown  in  Figure  4.6.  From  this  set  of  K  possible 


sampling  times,  any  selected  set  of  N  sampling  times  would  yield 
a  sampling  spectrum  of  period  equal  to  K/m,  where  m  is  some 
positive  integer.  We  are  interested  here  only  in  those  sets  of 
sampling  times  for  which  m  -  1.  These  will  be  the  sets  satisfying 
Equation  (4.15);  i.e.,  those  sets  for  which  K  is  equal  to  the 
smallest  common  denominator  K'  of  the  N  equations  given  by 
Equation  (4.15).  Of  course,  the  actual  values  of  the  spectral 
coefficients  a^  depend  on  the  specific  set  of  N  sampling  times 
chosen.  Since  the  number  of  possible  unique  sets  of  N  sample  times 
from  K  possible  sampling  times  (the  number  of  combinations  of  K 


things  taken  N  at  a  time  is 


is  n:  (kIn)  : 

w  — 


)  becomes  very  large  for 


large  K  and  N  with  K  »  N,  additional  criteria  beyond 
periodicity  may  be  necessary  to  be  able  to  select  an  optimum  set  of 
N  sampling  instants. 
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Figure  k  .6  Division  of  Block  Interval  T  into  K 
Possible  Sampling  Instants  p 
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Consider  any  nonuniform  block  stapling  function  slallar  to  that 
of  Figure  1».3  with  N  samples  per  block  Interval  T^,  and  with  a 
periodicity  K.  The  apectrua  of  each  such  sampling  function  Is 
completely  described  by  a  knowledge  of  the  c  xsplex  spectral  c  “efficients 
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Afaln,  t he  perl41clly  relation  of  Equation  (4.12) 

Intenchan^ln*  the  order  of  MMtl  >n,  Equation  (4.16' 


1 

1 


r. 


4  1 


1 

B 


•  J2«l 


n-1  4-1 


N 

*.n 

r2 

• 

1 

i  - 1  1 

.i  1 


2*H 


1  •  • 


4.19) 
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ih*  largest  and  smallest  spectral  components  is  noticeable.  As  will 
be  seen  later,  the  presence  of  excessively  large  spectral  co^jonents 
aajr  be  undesirable  in  some  applications  and  a  more  uni  fora  distribution 
of  power  versus  frequency  required. 

The  complex  conji^ate  sywetry  about  the  midpoint  k  -  10  given 
by  Equation  (4.14)  is  evident  in  all  the  exa^les  of  Figure  4.7.  The 


2 

.fed  values  jf  a^ 


as  shown  on  Figure  4.7  for  each  example,  can  be  compared  with  the 

theoretical  values  from  Equation  (4.17)  ana  (4.20)  of  =  jj  =  0.1 

for  all  examples  and  a^  0  for  Example  (d)  and  a£  =  =  0.1 

for  the  remaining  examplea  since  mr  =  l  for  these  examples.  All 

P 

computed  values  are  seen  to  be  within  ♦  O.ji  of  the  theoretical 
values--well  within  computational  accuracy. 

Another  important  special  case  of  nonunifom  sampling  (see 
Figure  4.6)  is  that  of  moving  every  second  sample  of  a  uniform 
sampling  function  by  a  constant  amount.  For  example,  this  could 


describe  the  quadrature  sampling  process  in  Chapter  II  and  III  of 
taking  quadrature  pairs  of  samples  at  a  uniform  rate.  This  example 
can  be  assumed  to  be  the  sum  >f  two  uniform  sampling  functions  with  a 
fixed  time  origin  offset.  Accordingly,  the  frequency  spectrum  for  this 
sampling  function  is  equal  to  the  vector  sum  of  the  contributions  from 
each  uniform  sampling  function.  For  the  example  of  Figure  4.6,  the 
uniform  spacing  between  corresponding  samples  is  0.2  Tp  which,  in 


general,  results  in  a  non-zero  component  for  k 
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-  multiples  of  5* 
The  zero  spectral  component  at  k  10  occurs  because  the  vector 
contributions  of  the  two  uniform  sampling  function  cancel  each  other. 
The  periodicity  of  K  =  20  follows  from  Equation  (4  .  15)  since  there 

are  20  possible  sampling  instants  per  block  and  N  *  10.  Performing 
the  following  sums  over  one  spectral  period: 


k  1 


2 


and 


T 


0 


0.1 


likewise  is  consistent  with  Equations  (4.17)  and  (4.20) . 

The  concept  of  the  complex  frequency  spectrum  a^,  as  a 
vector  sum  of  complex  spectral  contributions  from  two  uniform  sampling 
functions  in  the  example  of  Figure  4.8  can  be  extended  in  the  general 
case  of  Figures  4.3  and  4.4  to  the  sum  of  contributions  from  N 
uniform  sampling  functions.  This  is  equivalent  to  considering  a^ 
as  a  resultant  vector  equal  to  the  average  of  N  unit  vectors  with 
varying  phases.  The  maximur.  value  of  a^  -  1  occurs  when  all  N 

vectors  have  the  same  (zero)  phase.  If  the  N  vector  contributions 
are  uniformly  distributed  in  phase  over  2n  radians,  the  resultant 
a.  is  zero  [e.g.,  see  spectrum  of  Figure  4.7(a)  for  1  <  k  <  9). 
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4.6  Minimizing  the  Spectrum  Variance  by  Choice  of  Sample  Times 

In  Section  4.4,  it  was  noted  that  although  the  average  power 

2 

a^  of  the  sampling  spectrum  components  is  a  constant  which 

depends  only  upon  the  number  N  of  samples  per  block,  the  variance 
2 

of  a,  is  strongly  a  function  of  the  specific  sampling  times 

K  2 
In  this  section,  it  will  be  shown  that  the  variance  of  a^ 

can  be  minimized  by  proper  choice  of  £t  ^  and  the  conditions  on 

for  minimum  variance  will  be  derived. 

The  power  associated  with  each  a^  is  given  by 
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N  N 
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n^m 

2  1 

Since  the  mean  value  of  was  found  to  be  equal  to  ^  in 

2  2 

Equation  (4.17),  the  variance  of  a^  over  one  period  of  k 

can  be  defined  as: 
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^  1 

Substituting  the  expression  for  (  a^  -  ^)  from  Equation  (4.21) 


into  Equation  (4.22)  gives  as: 
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If  the  terms  for  which  (tn-tm)  =  (t^t . )  are  isolated  and  the 
summation  over  k  is  interchanged  with  the  other  four  summations, 
Equation  (4.23)  becomes: 


N 


N 

V 

/ 


2 

?  nti  „Tl 
n/m 


1 

K 


K 

\ 

i_ 

k=I 


cos 


N  N  N  N 


2  \  \  \ 

+  ~T  l.  L  i_  — 

N  n=i  m=l  i=l 


n/m 

i/j 

(i,j  )/(n,m) 
(i,j)/(m,n) 


(t  -t  )1 

T  v  n  m' 

P  J 


1 

K 


V 


cos 


k-1 


“•  <W 


(It.  21+) 


The  conditions  for  minimizing  by  choice  of  ^tn^  can  be  found 

by  examining  the  two  bracketed  sums  of  Equation  (4.24).  Let 


X 


I 

K 


cos 


2  2irk 


(t  -t  ) 
K  n  nr 


k=l 
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and 


_K 

1  \  2nk  /.  ,  v  2nk  N 

K  /  COST“  (V^*008—  (W' 

kTl  P  P 


Using  trigonometric  identities,  and  S  can  be  written  as: 


1  +  i- 

2  2K 


k=l 


4nk  /.  .  v 

cos  — —  (t  -t  ) 
T  '  n  m' 
P 


(4.25) 


and 


1_ 

2K 


T 

k=l 


2itk  ,  -  .  \ 

cos  ■= —  (t  -t  +  t. -t . ) 
T  '  n  m  i  J 
P 


+  4?  /  COS  (t  -t  -  t.+t.). 

2K  £_  T  v  n  m  l  j 

k-1  P 


(4.26) 


An  examination  of  Equations  (4.25)  and  (4.26)  shows  that  because  of 

the  periodicity  of  Equation  (4.12),  the  three  summations  over  k  are 

always  either  zero  or  positive  depending  upon  the  choices  of  tn,  tm 

2 

t. ,  and  t..  Therefore,  minimizing  o„  is  equivalent  to  minimizing 

1  j  K 

and  Sg  separately  over  all  sampling  instants. 

First,  consider  as  given  by  Equation  (4.25).  Using  the 

periodicity  condition  of  Equation  (4.12)  in  an  analogous  manner  to  that 
done  in  Section  4.4,  it  can  be  shown  that 


108 


1 

K 


_K 

\ 

/_ 

k=l 


hnk  /,  .  . 

cos  - —  (t  -t  ) 
T  n  m 

P 


{ 


1  whenever  (t  -t  )  -  + 

v  n  nr  —  2 

0  otherwise. 


(4.27) 


Equation  (It  - 27 )  shows  that  to  minimize  av  ,  the  sampling  times  should 
be  selected  to  minimize  the  number  of  occurrences  of  the  condition 


(Vtm)  =  2*  ’ 


(4.28) 


where  t  and  t  ai-e  any  two  sampling  times  with  t  >  t  . 
n  m  J  *  n  m 

Next,  consider  the  sum  given  by  Equation  (4.26).  Again, 

using  the  periodicity  relationship  of  Equation  (4.12),  the  two 
equations  analogous  to  Equation  (4.27)  are 


fl  whenever  (t  -t  +t.-t.)  =  -T  , 
1  v  n  m  l  j  p’ 


1 

K 


2nk 


k=l 


cos  ^r-  (t  -t  +  t  -t  )  ={ 

I  n  m  i  j  \ 

p 


0,  or  T  , 
P 


0  otherwise, 


(4.29) 


and 


K 

1  V 


fl  whenever  (t  -t  -t.-t.)  =  -T  , 
1  n  m  i  j '  p’ 


Znk 


K  /  C0S  Y-  ( +  Vt  )  =< 

fel  P 


0,  or  T  , 
P 


0  otherwise. 


(4.30) 


To  minimize  ,  the  sampling  times  should  be  selected  to  minimize 

the  number  of  occurrences  of  the  conditions  given  by  Equations  (4.29) 


and  (4.30).  These  two  conditions  on  t  ,  t  ,  t. ,  and  t.  can  be 

n’  m*  i’  j 
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r 


combined  into  an  alternate  equivalent  criterion  that  covers  all  cases 
and  is  easier  to  compute  in  general.  Consider  any  four  sampling 
instants  t^,  t^,  t^,  t^  ,  at  most,  two  of  which  can  be  equal,  order 
them  such  that 

tL  <  tg  <  t3  <  t1(  (J*.31) 

and  form  the  six  linear  combinations 

h  =  -  V  i  ^2  -  V> 

l2  -  - V  i  (t3  -  V 

and 

L3  =  "  ti)  i  ^3  ‘  V  •  (^-32) 

Then,  from  Equation  (h . 2 6),  Sg  is  increased  by  l/2  whenever 
L..  ,  L  ,  or  L,  equals  T  or  zero. 

L  u  p 

For  the  case  above  where  two  of  the  four  sampling  instants  in 
Equation  (U.31)  are  equal,  Equations  (U.31)  and  (U.32)  reduce  to: 

ti  <  tg  <  tj  (U  .33) 

and 

L1  '  (t3  -  t2)  1  (t2  '  V> 

l2  -  (t3  -  t2)  ±  (t3  t3) 

and 

L3  =  (t3  "  V  ±  (t2  "  tl)  •  (l4’3l+) 


2 

To  summarize,  in  order  to  minimize  ,  it  is  necessary  to 

minimize  the  number  of  occurrences  of  the  following  three  conditions: 


T 


(1) 

(t  -  t  )  =  7^  ,  for  t  >  t  , 

n  m'  2  ’  n  m’ 

(*i.35) 

(2) 

L,  ,  L0 ,  or  L,  r-  T 

1*  2’  3  p 

(14.36) 

(3) 

1*1,  Lg,  or  L3  =  0, 

(*.37) 

where  L^,  Lg,  are  given  by  Equations  (14.32)  or  (I1.3I4). 

Consider  a  given  sampling  function  ^t^  consisting  of  N 

distinct  sampling  instants.  In  order  to  determine  the  total  number  of 
occurrences  of  the  conditions  of  Equations  (4 . 35 )  -  (I4.37)  for  all 
combinations  of  sampling  instants,  the  following  procedure  can  be  used: 

(1)  Form  all  possible  pairs  of  two  distinct  sample  points 

t  and  t  such  that  t  >  t  and,  for  each  pair, 

compute  the  difference  (t  -  t  ).  Count  the  total 

number  of  occurrences  Mg  of  the  condition  of 

Equation  (U . 35 ) ;  i.e.,  the  number  of  times  that  (t  -  t  ) 

n  m 


(2)  Form  all  possible  sets  of  three  distinct  sample  points 
t^,  tg,  and  t^  and  order  each  set  according  to 
Equation  (14.33).  For  each  set,  compute  the  six  numbers 
represented  by  L^,  Lg,  and  of  Equation  (Ij.34).  Count 
the  total  number  of  occurrences  Mj  of  either 
Equation  (it .36 )  or  Equation  (1).37);  i.e.,  the  number  of 
times  that  ,  Lg,  or  equals  0  or  T  . 


Ill 


(3)  Form  all  possible  sets  of  four  distinct  sample  points  t^, 
tg,  tj,  t)(  and  order  each  set  according  to  Equation  ( . 3 1 ) . 
For  each  set,  compute  the  six  numbers  represented  by 
L^,  Lg,  and  of  Equation  (1.32).  Count  the  total 
number  of  occurrences  of  either  Equation  (1..36) 
or  Equation  (^ . 3T ) - 

If  no  occurrences  of  the  conditions  of  Equations  (I4.35)  -  (h .37 ) 

are  found  for  all  possible  combinations  of  sampling  times;  i.e.,  if 

Mg  =  Mj  =  M^  =  0,  then  si  ~  f  »  S2  0  and  the  minimum 

2 

value  of  the  variance  is  attained.  From  Equation  (I4.2I4)*  this 

value  is: 


/  2v  2_ 

°K  'min  N2 


(14.38) 


Equation  (h . 38 )  is  useful  as  a  standard  against  which  the  variance  of 
any  non-optimum  (i.e.,  non-minimum  variance)  sampling  scheme  can  be 
compared. 

2 

A  simplified  expression  for  c?K  from  Equation^  .2*4 )  for  any 
non-optimum  sampling  scheme  can  be  written  in  terms  of  Mg,  M^,  and 


M. 


An  examination  of  Equation  (h . 2h )  shows  that  each  Mg  violation 


results  in  two  non-zero  terms  and  each  Mj  or  M^  violation 

2 

results  in  four  non-zero  Sg  terms.  Thus,  for  the  general 

case  in  terms  of  Mg,  M^,  and  M^  is: 


(*.39) 


’k2  7  (1  •  7p  +  T  (V  4  T  (M3  +  V 

I.  h  I< 


(’K2)in  *5: 


The  examples  of  Figure  1*.9  illustrate  the  variance  minimization 
concepts  discussed  in  this  section.  Example  (l)  shows  one  of  many 
possible  sample  time  choices  (for  N  -  i*  samples  per  block  and  20 
possible  samling  times  per  block)  which  yield  a  minimum  variance 
spectrum.  Since  N  =  U  for  these  examples,  the  minimum  variance 


from  Equation  (U.38)  is  (c„  )  . 

K  rnn 


7  <l-5>  -  oM9- 


o 

.  ► 


Examples  (2)  -  (1*)  show  some  typical  non-minimum  variance  cases.  The 

table  at  the  bottom  of  Figure  *.9  gives  the  valu  s  of  Mg,  Mj,  and 

for  each  example  and  compares  the  computed  variance  with  the  theoretical 

variance  from  Equation  (U.39).  In  all  cases,  the  computed  variance 

was  accurate  to  within  computational  accuracy.  Note  that  the  minimum 

variance  spectrum  of  Example  (1),  while  not  perfectly  uniform  over  k, 

has  fewer  extremely  large  or  extremely  small  values  compared  with  the 

other  examples.  The  uniform  sampling  case,  Example  (5),  illustrates 

the  upper  extreme  of  maximum  variance  because  all  of  the  spectral 

power  is  concentrated  at  the  minimum  number  of  frequencies.  The  value 

of  oK  for  the  uniform  sampling  case  is  given  by  oK  =  l/N  (1  -  l/N). 

2 

Therefore,  the  range  of  is  given  by: 
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EXAMPLE 

M2 

m3 

m4 

COMPETED 

THEORETICAL 

o-K* 

1 

0 

0 

0 

04665 

0469  (MIN.  VARIANCE) 

2 

1 

1 

0 

0705 

0703 

3 

1 

0 

2 

0867 

0859 

4 

1 

2 

0 

0860 

0859 

5 

2 

4 

4 

.1875 

.1873  (MAX.  VARIANCE) 

Figure  4.9  Examples  of  Variance  Minimization  Principles 


11U 


ft 


V 


I 


4  * 

«  ► 


) 


h  (1  *  ^  - 

N 


The  minimum  -ariance  sampling  points  shown  in  Example  (1)  of 
Figure  b  .9  were  relatively  easy  to  find  by  trial  and  error  because  of 
the  relatively  small  number  of  possible  sets  (eleven)  of  2,  3,  and  4 
samples  which  had  to  be  checked  against  the  conditions  of  Equations 
(4.35)  -  (4.37).  However,  as  the  number  of  samples  per  block  N 
becomes  large,  this  method  rapidly  becomes  unmanageable  except  perhaps 
with  a  digital  computer.  If  minimum  variance  solutions  do  exist  for 
large  N,  it  seems  probable  that  a  computer  algorithm  could  be  devised 
that  would  find  them.  Investigation  of  this  topic  is  left  as  a  possible 
area  for  further  research. 

Even  if  minimum  variance  solutions  cannot  readily  be  found  or 

are  not  practical  for  a  given  application,  the  minimum  variance 

2 

concept  is  important  because  the  minimum  oK  given  by  Equation  (4.38) 

2 

and  the  oK  for  the  uniform  sampling  case  give  two  extremes  against 
which  any  practical  sampling  scheme  can  be  compared.  This  will  be  done 
in  the  next  section  for  a  unique  nonuni  form  sampling  technique  which 
uses  a  pseudo-random  number  generator  to  select  the  sampling  times 


4.7  Sampling  Time  Selection  Using  Pseudo-Random  Number  Generators 
This  section  describes  a  nonuniform  sampling  scheme  which  uses 
a  pseudo-random  number  generator  [19]  to  select  the  nonuniform  sampling 


N 

times  ft .  j  .  Although  the  scheme  does  not  generate  a  minimum 
i  - 1 


! 
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variance  spampling  spectrun,  it  achieves  an  important  unique  spectrum 
shaping  by  suppressing  the  amplitude  of  the  ,,  spectral  components 
for  small  values  of  k,  it  is  easily  reproducible,  and  it  cam  be 
readily  modified  to  obtain  any  desired  number  of  samples  per  block, 

N.  In  this  section,  the  basic  pseudo-random  sampling  function  will 
be  analyzed. 

Let  each  block  sampli*'.  time  be  divided  into  N  equal 

segments  of  duration  T  =  T  /N  as  shown  in  Figure  ^ .  10(a).  One 

s  p 

sampling  instant  is  to  be  selected  at  some  point  in  each  T  interval 

s 

to  give  the  required  total  of  N  sampling  times  per  block  interval. 

Now,  let  each  segment  of  length  T  be  subdivided  into  N  uniformly 

s 

spaced  possible  sampling  instants  spaced  A  =  T  /N  apart  as  shown 

s 

in  Figure  U . 10(b).  Only  one  of  these  N  possible  sampling  instants 

will  be  selected  by  the  sampling  logic  during  each  T  second  time 

s 

segment.  The  choice  of  which  of  the  N  possible  sampling  instants  to 

select  during  each  T  segment  is  made  by  a  pseudo-random  number 

s 

generator  as  shown  in  Figure  In  this  figure,  a  pulse  train  with 

a  period  of  A  seconds  between  pulses  drives  a  binary  counter  which 

counts  to  N  in  Tg  =  NA  seconds  before  being  reset  to  zero  to 

repeat  the  count  during  each  Tg  second  interval.  The  Pseudo-Random 

Number  Generator  (FRNG)  which  is  assumed  to  have  a  repetition  cycle 

of  N  states  is  driven  by  the  pulse  train  divided  by  N.  The  PRNG 

thus  changes  states  at  the  beginning  of  each  T  second  segment  shown 

in  Figure  U.  10(b)  at  the  same  time  at  which  the  binary  counter  is  reset. 

At  some  time  during  each  T  second  interval  in  which  the  PRNC  remains 
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BLOCK  INTERVAL  *  Tp  *  NT, 


. ♦  — +-  ♦  «  4  "•  t  — - 4  — ♦ - f 

0  T,  2T,  3T,  •  •  •  nT,  —  Tp*NT,  t 

(o)  SUBDIVISION  OF  BLOCK  INTERVAL.  T>t 
INTO  N  EQUAL  SEGMENTS 


T,  •  NA 


III! - — - >  * - 1 - 1  - 

0  A  2A  3A  •  •  •  nA  •  •  •  T,»NA  » 

( b )  SUBDIVISION  OF  EACH  SEGMENT  OF  LENGTH  T, 
INTO  N  POSSIBLE  SAMPLING  TIMES 


Figure  I*.  10  Construction  of  Possible  Sampling  Times  for 
Pseudo-Random  Sampling 
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N 


f  (vo^niforw  sampling  times  ( t,  }  is  generated  In 

^  J  i-l 

•  ft  t  seconds  (one  block  Inteval)  and  the  entire  process  is 
repeated  fndeflnltely  since  the  TONG  is  periodic. 

ft*  HOW  could  be  connected  in  a  variety  of  ways  to  give 
sequences  of  varying  lengths.  For  example,  the  PRNG  feedback  could  be 
connected  to  give  a  maximal  length  sequence  for  which  the  theory  has 
been  well  developed  (19).  If  the  PRNG  has  an  m  stage  shift  register, 
the  maximal  length  code  length  is  (2m  -  l).  Thus,  by  making  the 
block  length  N  equal  to  the  code  length  (2m  -  l),  an  effective 
sampling  time  generator  can  be  made  whose  sampling  spectrum,  although 
not  of  minimum  variance,  exhibits  a  high  degree  of  randomness.  With 
minor  modifications,  a  PRNG  can  be  made  to  generate  shorter  repetitive 
sequences  of  varying  lengths,  thus  giving  more  freedom  in  the  selection 
of  a  sampling  scheme. 

Assuming  the  system  shown  in  Figure  4.11  is  used  to  generate  a 
set  of  repetitive  nonuniform  sampling  times  with  N  samples  per 
block,  the  results  of  this  chapter  can  be  used  to  predict  the  properties 
of  the  sampling  spectrum.  The  expression  for  the  sampling  coefficients 
a^  can  be  found  by  letting 
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in  this  case ,  sin.  <?  ?  T  /K  and  A  T  /H  .  Jlnc*  ih*>  t>loek 

5  p  •’ 

Interval  T  is  subdivided  into  N2  possible  equally  spaced  stapling 
times  as  shown  in  Figure  1.10,  the  t-vrl  dlcity  K  >f  as  given 

by  Equation  (1 .  _c' '•  i"  K  -  N  and  the  spectral  period  in  Ht  is 

2  i. 

N  f  ,  where  f  =  .  The  average  power  of  a^  as  given  by 

P  P  P 

Equation  (1.17)  is  1.  constant  equal  to 


1 

K 


1 

N 


and  the  average  value  of  as  given  in  Equation  (1.20)  is  either 

l/N  or  0  depending  on  whether  or  not  the  last  sample  t^  occurs 

at  the  end  point  of  the  block  interval.  T  . 

P 

As  an  example  of  this  technique,  a  computer  program  was  written 

no  simulate  the  sampling  time  generator  of  Figure  1.11.  A  FRNG  with 

a  J.ut&ge  shift  register  was  simulated  with  feedback  connections  to 

3 

give  a  maximal  length  sequence  of  length  N  =  7  =  2  -  1. 

Figure  1.12(a)  shows  that  each  basic  block  interval  of  length 
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fttfwM  (MW  I  ■#»!  (♦  ;*•  •***  fti-*ft  *■»-(•* ‘*4  ,  (ft*  *'<•«»>*(*(  I4»-  of  ^ 

«»IWU(«  of  ft***ftdt*d  ft  *•*(«•♦  ♦»*«*  •***•  *ft»(o*  (*  (. 

ft  *****  *«*»•.*  “f  *t#(a  ft**  ♦**;;  ft,  <44  •  f  »<•««*  '*•  * 

>«r(**»  <-  ft*  fft(*W  |UdUMM  “'»**  (*•  fft  t*ft*« 
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(ft*  •***,*•  *M4  *  *o«*k«*(«  **><**-<  of  •’«$$#*>••  i oft.  i*  *Mc<*4  wh*ft. 

*  ^  ft,  «*•(*  ***(.«#  ift»  *  f>;;  lift*  *****  of  }♦,  (t* 

*ff**(  U  <mo  toft***  »«••**-(  ft*4  ft*ft?t*Hft  fl  «*(  ,*<  l  o*„  of  ft^  !*fM  * 

U  o***r»*#„  "U*  «ff*r(  •lit  ft*  »*(■.**  |ttft*(rft(*4  }ft*.*r 

Ifc  fti#***  ft  1*  *****  *  l*oni  •aaftil**  *j>**,r*»  U  j>t»((*4  for  s  -  *>„ 

it.  CH«j(*f  i,  (ft;*  ftffft<(  till  ft*  \>  fftt*  iftjwriwi  lafiieftUMi 

I*  (»*  of  44  tftfftoftU  f*ip3fvi*i  for  correlation  f^U^t 

0 *t  c,lpf>*4  tlrv*ftoi4ft; 

7>*  ***at>l la*  fft*x(|im  Ihoftf.  in  figure  ft.U'c)  for  Care  6,  which 
tfti  j*.\f  j*w  jccvrr*«*  of  **ro  in  Equation  ft.J6)  or  Equation  (ft, 37), 
*%*  (h*  MX  *.h*(  c j'jld  b*  found  by  •  Halted  trial  and  error  starch. 
Th*  question  of  existence  of  a  tra#  nlnlau*  varlanc*  ■  sapling  function 
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ft1)  Is  still  op«n 
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<A4MWC<C<UIM  f4o«4CU«f  »y  a; I  HWti  teftc«A<e*  HI|U*i  At'* 

*  «^9  |  *  4  tAC  IAU1MM  »f  <fCC<<%|  <i(W|l.ACAA»  ** 

A  -  -I  *  “j  U  <*U  «♦*•»-  51MA  *<**»<*<  M>-W  f'1"* 

*  <;  t\|  ;e*U-*C  A  A  *,  «*eie  CUM-  M<e*c* 

#»*»** 


f«<  *;;  «  tw*4M  *  «***•*  a;;  u<«*e<  *%;•**  u  <*e 

MUf»A;  ;  >  *.  -  *  ^  a^wnaMoa  «f  fcf«c<io*.  a  a#  <aa  $• 

CoA4l4c<MI  l.  M  <AC  AMIU4  of  A  rtfUf*  Aietfitoted 

M  f***c  o»*f  }»  AMI  M  '*>«*  M4l  <»  •*<» 

A  CJAfcrUoA  of  CAACA  l -a  IkMA  U4I  <*C  *J4WT<r«  1  •»*{*  «4 
nrlAC<  t*  wy  *»«.  "Me  efcolM  of  MliMt  >Wi  «'4». 

!>»4,  Ifc  A  prACliCAl  frfclee,  All  (DM* IAEA  rode#  »I»w44  Ac  *HAln<4 
in  jrder  to  etke  a  ~beet  choice.  The  tAffpling  code  for  Cam  *>  it 
ldenticel  to  thet  of  Caaa  i  escept  thet  a  different  initiAl  iterting 
point  in  th«  sequence  mi  tssuMod,  which  it  equivalent  to  t 
redefinition  of  the  tin*  origin  of  the  stapling  block.  The  vtrlence 
And  All  apcctrtl  components  |t^|  were  found  to  be  identical 
Independent  of  the  starting  point. 

The  results  of  the  cospArlton  of  the  exeaplet  of  Figure  1*.12 


ctn  be  aunntrlxed  At  follows: 
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*  thMi  4m«(VH  »iwu*  «4*l*H  ♦  #•**;;  *  *4f;n.*4* 
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(i  *  4-  ««*«***•*  i+h  (.«  af  ♦;;  iim»  4*»««44«#  ♦$•**#<* 

«M*  *  %*  U  u*1  •;;  **ff«4«4l  *<**♦„***.<,♦  ** 

4  *  ««4*|»$S*  af  *  4*« 

«  *»*..  fhttt^ifar******  *444 ;  1 44  F,**4|o**  <«t  *•  f'<*«*M#  *»« 

;^h»#  .  *•♦*  *1*4*4*  *4*  <44k«  **ffU*4  !•♦.*«*  |  »**, 

ir  i*<r/  «*u*.  *i«  f#  rt*4 

******  *«4*tl4  »f  4  fW|Utf  r,*<4|*»*- 

U  */  f}#+#«  4»|i  *M«*  f*»a*4  (t-MUl  **°*M o*t  of  *** 

•4Mflt*4  *p+* 1  rjf  4»t«r*'*4  4;/-  4  *,4*4tf»  Mwxi  t%3 Hi  c>m4*  af 

Uftftfc  t>  -  l>  *  If*  U  -«M4  l»  Mlwf  «4|(*  Of  If*  |M**l*i4 

««*pllf>*  11 14*4  U  444*  ?#  |M4T#4l  M*  fl#+»  4.|0t»f!  14  Cfca*#*.,, 
Mo«9V4r.  Jftljf  s  -  ii  Mffltl  f*f  *l«l  |M#r#4l  *f*  4*4,  40  'M 

_  a  j|» 

HOG  If  r*M(  ifitr  ’.I*  *  *  njMr  Ia  it*  co4»  *4  ft*  »«ftl«4 
procMf  If  r4p44t4d.  Th#  ifDffitnf  *p*e* rm  ’.l»4  mi  4  p»rlo4 
k  -  /’ill?)  II,  1 5*  Which  If  #9*41  l«  IM  VoUl  As*4#F  of 

pofflbl#  fiMplinf  Umm.  Although  only  4CojI  1/12  of  4  ion:  ip*ci  r*l 
ptrlod  If  fhown  In  Flfur*  *.  13,  th*  ikUo <  fra*  k  - 


IOO  vo 


k  -  &W  1#  i/pieki  of  iK#  general  t>#|*  >f  <*#  +isfl<>u*4 

portion  of  ihe  ipwinn.  l!iw#  #aeh  spectral  component  ^  u  f»r«N 
by  averaging  a  large  waktr<f  unit  vtcun  wilh  HiMUiUjf  rv4a* 
phuM,  on#  would  nupwt  Umi(  If  ihe  InltUl  regl<*»  for  I  ^  k  <  » 
la  eliminated,  ih#  probability  distribution  function  for  u#  reaeiMf* 
•pootral  components  |a^  |  should  b#  approximately  haylelgh  with  in# 
following  par— atara  for  h  «  96: 

*  i  «  O.OIOH7 


»  *  0.0905. 


This  waa  experimentally  rhovn  tc  :*  a  valid  assumption  for  larga  M. 

Tha  following  tabl#  cowparaa  ttr  r«asur*d  spectrua  parameters  shown 
for  this  example  with  tha  theor**lcal  valuas  darlvad  In  this  chapter. 


Parameter 


V.-asured 

Value 


Theoretical 

Value 


K 

1  v" 
K  t 
k-1 


0  010J1182 


jjj  =  O.OIOI1I67 


O.rr  01818 


'“/W  0-0°°W 

(a2) _ =  0.010308 

f\  alclX 


!*•  preffl-t  of  I^f9«*|wv  for  MMll  fc  il  Mil 

lll»*tr*t*4  Vf  of  figure  k ,  IS  which  chow*  on  eppruslaete 

I  In**/  growth  of  jw^J  for  l  ^  li  <  96  with  ih*  region  of  the 
«*st*4»  suppression  •ffrcUmtMi  for  l  *  k  <  kB.  The  ieporlsnce  of 
in*  mm l  i  k  sup* cession  effect  In  the  reduction  of  odd  hemonlc 
r«*t>>Ai«i  In  e-rreletijn  functions  of  clipped  sinusoldel  slgnels  trill 
t>e  discussed  In  the  nest  chepter. 


CHAPTER  V 


APPLICATION  OF  BLOCK  SAMPLING  TO  CORRELATION 
FUNCTIONS  OF  CLIPPED  SINUSOIDAL  SIGNALS 

5.1  Introduction 

In  Chapter  III,  correlation  function  properties  of  clipped 

sinusoidal  signals  with  uniform  sampling  were  investigated.  As  shown 

in  Figure  3.10,  the  uniform  sampling  produced  undesirably  large 

correlations  for  some  reference  frequencies  due  to  the  presence  of 

odd  harmonics  caused  by  clipping.  For  example,  the  frequencies  that 

respond  to  the  third  harmonics  have  nominal  correlations  of  l/3 

of  the  fundamental  correlation.  Likewise,  the  fifth  harmonics 

generate  correlations  of  about  1/5  of  the  fundamental  response  at 

s  t 

some  frequencies  and,  in  general,  the  (2n  +  1)  harmonic  produces  a 
correlation  of  about  ^  relative  to  the  fundamental  response. 

In  Section  3.5,  the  use  of  nonuniform  "block"  sampling  was  proposed 
to  reduce  these  odd  harmonic  correlations  by  destroying  the  synchronism 
between  the  uniform  sampling  peridocity  and  the  periodicity  of  the  odd 
signal  harmonics. 

In  Chapter  IV,  'block  sampling"  was  defined  and  some  of  its 
properties  were  developed;  namely,  the  spectral  periodicity  properties, 
the  average  power  spectrum  properties,  the  requirements  for  minimum 
variance  of  the  sampling  coefficient  power  spectrum  and  small  k 
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spectral  amplitude  suppression  for  FRNG  sampling.  In  this  chapter, 
these  block  sampling  results  will  be  applied  to  the  reduction  of  odd 
harmonic  responses  in  correlation  functions  of  clipped  sinusoids.  A 
specific  correlation  function  example  is  included  using  a  FRNG 
generated  sampling  function  similar  to  the  one  discussed  in  Section  U.7. 

This  chapter  considers  the  basic  reasons  for  the  ability  of 
block  sampling  to  reduce  the  odd  harmonic  responses  of  correlation 
functions  of  clipped  sinusoidal  signals.  In  the  next  section,  a 
qualitative  correlation  function  analysis  with  block  sampling  will  be 
made  to  help  develop  insight  into  the  problem  and  to  help  predict 
which  of  the  reference  frequencies  will  have  the  greatest  response 
for  a  given  input  frequency. 


5.2  Qualitative  Analysis  of  Correlation  Functions  with  Block 
Sampling 

The  ideal  block  sampling  quadrature  correlation  function  model 

of  Figure  5.1  is  useful  for  analyzing  the  relationships  between  the 

input  frequency  f^  and  the  reference  frequencies  fR  which  produce 

significant  correlations.  Since  the  frequency  difference  (fR  -  fj) 

is  the  important  parameter  in  this  chapter,  we  will  adopt  the  point 

of  view  of  holding  the  input  frequency  f^  constant  while  varying 

the  reference  frequency  fR  to  study  the  correlation  function 

dependence  upon  the  frequencies  fj  and  fR.  The  spectra  of  the 

pre-envelopes  p  (t)  and  p  *(t-T)  at  the  sampler  inputs  are  shown 
x  y 

in  Figure  3.6,  while  the  spectrum  of  a  typical  block  sampling  function 


<*■ 
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is  shown  in  Figure  1* .5.  The  spectrums  of  the  sampler  outputs  p  (t^) 
and  p  *(t .  -t  )  are  equal  to  the  convolution  of  the  sampler  input 

y  * 

spectrum  (Figure  3.6)  and  the  sampling  spectrum  (Figure  4.5).  Thus, 
if  X(f)  and  Y*(f)  represent  the  Fourier  Transforms  of  px(t)  and 
p  *(t)  and  P(f)  is  the  block  sampling  spectrum  given  by 

y 


Equation  (4.7) >  then 


oo 

x(f)  =  1  T  T^-TT  p[f  -  <-l)n  (2n  +  *>  V 


n=0 


kfp], 


n=0  k=-o 


(5.1) 


OO 

Y*(f)  '  7  T  f[f  M-Dq  (2q+ 1)  fH] 


q=0 


OO  or 


■III 


(-l)q  a£ 

(2q  +  1) 


q=0  £=-» 


Btf  +  (-1)Q  (2q+  1)  fR  -  It  ], 


(5-2) 


where 


a^  =  sampling  spectrum  coefficient  (see  Equation  1*.6),  and 
Tp  =  l/fp  =  block  sampling  period  in  seconds. 


The  complex  spectrum  Z^(f)  of  the  correlation  function  z^(x) 
of  Figure  5.1  is  given  by: 


x(f)  Y-(-f) 
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k  -®  l 3 -m 

•  6[f  -  (-1)° 

(2n  +  1)  fj 

-  -Ip)  Mf  -  (-l)q  (2q  • 

»>  Vr ' 

P 

V 

'5.J ) 

Although  Equation  (5.3)  Is  not  a  useful  analytic  expression  for 
computing  the  actual  correlation  function,  It  Is  useful  for  determining 
which  reference  frequencies  fg,  combined  with  a  given  Input  frequenqr 
fj.,  can  cause  the  largest  odd  harmonic  correlations. 

Typical  plots  of  X(f)  and  Y*(-f)  from  Equations  (5.1)  and 

(5.2)  are  shown  in  Figures  5*2  and  5.3.  These  plots  show  that  the 

total  spectrum  of  each  sampler  output  consists  of  the  sum  of 

contributions  from  all  odd  harmonics  (only  the  fundamental,  third, 

and  fifth  of  which  are  shown  in  each  figure).  The  spectral 

t/h 

contribution  of  each  n  odd  harmonic  is  seen  from  these  figures 
to  be  the  complete  block  sampling  spectrum  of  Figure  4.5,  modified 
in  amplitude  by  and  centered  about  f  =  (-l)n  (2n».)  fj 

for  the  input  frequency  harmonics  and  about  f  -  (-1)^  (2q  *  1 )  fu 
for  the  reference  frequency  harmonics.  Each  spectral  component  can 
be  computed  using  Equation  (4.6)  for  a  specific  block  sampling 
function.  Thus,  all  of  the  properties  derived  for  block  sampling 


Fleur*  %3  Typical  Specl of  Saapled  Reference  Pre-Envelope 
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spectra  In  Chapter  IV  can  be  applied  to  the  analysis  of  correlation 
functions  for  clipped  sinusoidal  signals  with  block  sampling. 

Since  the  correlation  function  spectrum,  by  Equation  (5.3), 
is  equal  to  the  product  of  the  two  discrete  spectrums  shown  in 
Figures  5*2  and  5*3,  it  is  obvious  that  some  values  of  fT  and  fD 
will  produce  greatly  different  output  spectra  than  other  values. 

For  example,  values  of  fj  and  fR  can  be  picked  such  that  no 
frequency  components  of  X(f)  and  Y*(-f)  align,  and  the  resulting 
correlation  function  spectrum  (and  also  the  correlation  amplitude) 
will  be  zero.  On  the  other  hand,  for  some  combinations  of  f^,  fR, 
and  fp,  maximum  alignment  of  input  and  reference  spectral  components 
will  occur  which  will  tend  to  maximize  the  correlation.  More 
precisely,  since  each  component  of  the  block  sampling  spectrum 

is  a  complex  quantity,  the  spectrum  Z^(f)  consists  of  complex 
components  which  are  products  of  two  or  more  complex  coefficients. 
Since  the  correlation  function  is  proportional  to  the  integral 

OO 

I  Z_(f)  e^2n**T  df,  the  peak  correlation  depends  upon  the 

J-co  5 

relative  phases  and  amplitudes  of  all  of  the  complex  components  in 

Z^(f).  Although  precise  equations  will  t>e  derived  in  the  following 

sections,  it  is,  in  general,  qualitatively  true  that  values  of  f^, 

fD,  and  f  that  minimize  alignment  of  spectral  components  result  in 
n  p 

lower  average  correlations  than  values  which  maximize  spectral 
component  alignment. 

Using  this  as  a  rough  criterion,  we  can  derive  some 
relationships  between  fj,  fR,  and  f  which  result  in  maximum 
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alignment  of  spectral  components  and,  hence,  larger  average 
correlations.  Considering  first  the  spectrum  of  the  pre-envelope 
of  the  input  signal  as  shown  in  Figure  5.2,  it  can  be  seen  that  the 
relation 

^fI  =  mj  fp,  (5-M 

where  m^.  is  any  positive  integer,  is  a  sufficient  condition  to  insure 
alignment  of  all  input  odd  harmonic  spectral  components  on  a  set  of 
frequencies  with  f  separation.  For  this  special  case,  the  spectrum 
of  X(f)  reduces  to  one  similar  to  the  example  shown  in  Figure  5.4. 

In  the  same  manner,  the  relation 

ws  =  "eV  (5-5) 

where  nip  is  any  positive  integer,  assures  the  alignment  of  all 
reference  odd  harmonics  on  the  set  of  frequencies  shown  in  Figure  5.5. 
Each  of  the  complex  spectral  components  represented  in  Figures  5.4 
and  5.5  is  thus  the  vector  sum  of  an  infinite  number  of  weighted  complex 
coefficients,  one  for  each  odd  harmonic.  These  cases,  represented  by 
Equation  (5-4)  and  (5.5),  can  be  though  of  as  "Worst  Cases"  in  the 
sense  that  the  spectral  energy  is  distributed  over  the  fewest  possible 
frequencies  rather  than  being  distributed  over  a  greater  number  of 
frequencies  with  smaller  average  amplitude. 

Even  though  4fT  and  4f_  are  multiples  of  t  as  given  by 

A  n  p 

Equations  (5.4)  and  (5.5),  the  product  spectrum  will  be  zero  unless 


figure  5.4  Typical  Input  Pre-Envelope  Spectrum  When  4f. 


Figure  5.5  Typical  Reference  Pre-Envelope  Spectrum  When  4f. 
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the  non-zero  flpectve'  rcmponent*  of  fj  sh.<vn  In  Figure  **>..  align 
with  those  of  shown  in  figure  5.5;  i.e.,  unless 


(5.6) 


where  n.  <s  any  positive  or  negative  integer. 

Equations  (5.U),  (5. 5)  and  (5.6)  thus  define  the  conditions  on 
f . ,  fu,  and  f  which  result  in  the  Worst  Case  correlation  responses. 

in  p 

Let  the  input  frequency  f^  be  fixed  and  satisfy  Equation  (5.^), 
and  consider  the  problem  of  finding  all  the  reference  frequencies  fR, 
satisfying  Equations  (5*5)  and  (5*6). 

It  should  be  noted  chat  c1ust  because  a  reference  frequency 
satisfies  the  Worst  Case  criteria  docs  not  necessarily  imply  that  the 
corresponding  peak  correlation  wi.M  be  high.  The  actual  correlation 
amplitude  for  given  values  of  fj  and  fR  is  proportional  to  the 
vector  sum  of  a  large  number  of  complex  numbers  (spectral  coefficients) 
and  thus  can  be  relatively  small  or  relatively  large.  It  will  be  shown 
later,  however,  that  the  average  correlation  is  the  largest  for  the 
Worst  Case  frequency  conditions,  and  that  the  correlation  variance 
about  the  average  value  over  all  Worst  Case  frequencies  can  be 
relatively  large. 

If  fj  is  a  fixed  input  frequency  satisfying  i*f^  -  1*^, 

then  the  Worst  Case  reference  frequencies  satisfying  Equations  (5.5) 
and  (5-6)  are  given  by 


fn  =  i  i  ■  V 


(5.7) 
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where  m  1,  2,  3,  ...  .  Thus,  the  Worst  Case  reference  frequencies 
are  f  Hz  apart  for  a  fixed  input  frequency  -  j..  Each  time  the  input 
frequency  f^  is  'ncremented  by  \/k  f  from  i  ;s  initial  value,  a 
new  set  of  reference  frequencies,  each  separated  oy  f  Hz,  will 
satisfy  the  Worst  Case  conditions. 

It  i3  also  of  interest  to  briefly  consider  >ther  combinations 
of  fT,  f0  and  f  which  lead  to  a  non- zero  correlation  function 

1  X\  p 

spectrum  Z(f),  but  which  do  not  meet  the  Worst  Cas '  criteria.  For 
example,  we  can  let  the  input  frequency  satisfy  tlv.*  Worst  Case 
condition,  4f^  =  f  ,  and  ask  how  many  non-zerc  reference 

frequencies  exist  between  each  pair  of  Worst  Case  reference  frequencies 
given  by  Equation  (5-7)-  If  fj-  Is  held  fixed  and  fj.  is  moved  from 
its  Worst  Case  value  by  l/3  f^,  the  third  harmonic  spectral 
components  of  fR  (see  Figure  5.3)  will  align  with  the  ^'ctral 
components  of  f^  (see  Figure  5.U)  and  a  non-zero  correlation 
function  spectrum  will  result  whose  average  power  (as  wil  i  be  shown 
in  Section  5»*0  is  about  l/3  of  the  Worst  Case  average  power.  Moving 
fp  by  another  increment  of  l/3  f^  produces  a  similar  result.  Using 
a  similar  argument  for  the  higher  order  odd  harmonics,  the  set  of 
non-zero  response  frequencies  shown  in  Figure  5.6  can  be  con  tructed. 
Although  only  the  non-zero  reference  frequencies  caused  by  the  third, 
fifth,  and  seventh  harmonics  are  shown  in  Figure  5 similer  responses 
due  to  all  higher  order  odd  harmonics  are  also  present. 

Now,  consider  the  case  where  the  input  frequency  does  pjt 
satisfy  the  Worst  Case  condition  4fj  -  f  .  Let  f ^  be  ar 
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AVG.  CORR. 

RELATIVE  TO  NOTE:  ONLY  1.3, 5.7  HARMONIC  RESPONSES  SHOWN. 


fl  +  mfp  f|  +  (m+l)fp 


Figure  5.6  Non-Zero  Reference  Frequencies  When 
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arbitrary  distance  x  away  from  the  nearest  Worst  Case  frequency  as 
shown  in  Figure  5 -7(a)-  For  this  general  case,  the  input  and 
reference  harmonic  spectra  will  not  align  as  shown  in  Figures  5.4  and 
5.5,  but  will  be  similar  to  those  shown  in  Figures  5*2  and  5-3.  The 
odd  harmonic  reference  response  pattern  similar  to  that  of  Figure  5.6 
now  becomes  more  complex  because  the  alignment  of  all  possible  pairs 
of  reference  and  input  odd  harmonics  must  now  be  considered.  However, 
it  is  instructive  to  consider  the  most  significant  alignment 
combinations  as  shown  in  Figure  5.7(b).  This  figure  includes  all 
possible  alignment  combinations  of  the  fundamental,  third  and  fifth 
harmonics  of  the  input  and  reference  frequencies.  Figure  5 .7(b)  shows 
that,  relative  to  Figure  5.6,  two  more  reference  frequencies  respond 
to  both  the  third  and  fifth  harmonic  terms  and  that  the  average 
responses  of  Figure  5.7(b)  are  somewhat  lower  (by  about  20  percent) 
ohan  the  corresponding  responses  of  Figure  5.6. 

Figures  5.6  and  5.7  both  show  that,  although  minor  responses 
are  present  due  to  the  odd  harmonics,  the  dominant  terms  are  those 
for  which 

fH  *1  f  m  fp’ 

where  m  .i  l»  .1  2»  2  3»  ...  .  Although  the  Worst  Case  correlation 
responses  shown  in  Figure  5.6  are  somewhat  higher  (about  20  percent) 
on  the  average  than  the  corresponding  terms  in  Figure  5.7,  the 
correlation  variance  (fluctuation)  at  these  frequencies  will  be  shown 
to  be  large  enough  to  make  this  difference  less  significant. 


ll*3 


Careful  selection  of  the  sampling  spectrum  periodicity  K  is 
important  if  undesirably  large  extraneous  correlations,  equal  to  those 
for  the  uniform  sampling  case  of  Chapter  III,  are  to  be  avoided.  For 
example,  if  the  Worst  Case  conditions  are  met  and  if  the  third  harmonic 
periodicity  (see  Figures  5*2  and  5.3)  of  the  input  spectrum  coincides 
with  some  reference  frequency  fR,  then  that  reference  frequency  will 
produce  a  peak  correlation  of  about  33  percent  relative  to  the 
fundamental  response  of  unity.  That  is,  if  (-3fT  +  NKf  )  =  fD  for 

X.  P  A 

some  reference  frequency  fR,  where  N  is  any  integer,  then  this 
undesirable  situation  exists.  Similar  conditions  for  the  fifth 
harmonic  results  in  approximately  a  20  percent  correlation.  Since 
the  periodicity  K  is  usually  controllable  by  choice  of  sampling 
times,  undesirable  fundamental,  third,  and  fifth  harmonic  responses 
can  be  avoided  by  choosing  K  to  satisfy  the  following  inequalities: 

V""p  *  V 

fI  -  NK  fp  * 

fR  -  HK  fp  j  -Jfj, 

V^'p  * 

and 

VNKV  5fi 

for  any  f ^  or  fR  in  the  range  of  interest  where  N  is  any  positive 
or  negative  integer.  Similar  conditions  can  be  derived  for  higher 


harmonics . 
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The  qualitative  remarks  based  on  the  model  of  Figure  5-4  and 
the  frequency  plots  of  Figures  5*2  though  5-7  have  been  based  on  the 
assumption  of  sinusoidal  signals  of  infinite  time  duration  and  thus  of 
spectral  lines  of  zero  width.  For  signals  of  finite  duration  T, 
each  of  these  spectral  lines  have  a  finite  width  of  approximately 
l/T  Hz  and,  therefore,  the  spectral  picture  becomes  a  continuous 
"smeared"  one  instead  of  a  discrete  one  as  described  by  this  model. 

The  real  value  of  this  model  lies  in  its  ability  to  help  us  gain 
qualitative  insight  into  some  of  the  complexities  of  the  system 
response  and  to  help  predict  the  Worst  Case  reference  frequencies  for 
a  given  input  frequency. 

5.3  Deterministic  Analysis 

In  order  to  obtain  a  useful  quantitative  correlation  function 

expression  with  block  sampling,  consider  the  model  of  Figure  5.8. 

This  model,  based  on  a  similar  model  in  Tou  [36],  represent  nonuniform 

block  sampling  as  parallel  uniform  sampling  (at  the  block  rate  T  ) 

of  N  delayed  versions  of  the  signal,  where  the  delays  are  equal  to 

r  >N 

the  N  nonuniform  sampling  times  <  t.  >  in  each  block  interval. 

i  1 

w(mTp)  represents  the  average  of  the  N  samples  in  each  block 
interval  and  the  correlation  function  Zg(t)  is  formed  by  averaging 
over  a  signal  duration  time  of  M  blocks  or  T  seconds  where 
M  =  T/T  .  Since  the  input  and  reference  are  both  deterministic 
sinusoidal  signals  and  since  the  sampling  instants  ^t^  are  assumed 
known,  the  complex  output  correlation  Zg(i)  is  a  deterministic 
function  whose  real  and  imaginary  parts  are  equal  to  the  quadrature 


Quadrature  Correlation  Function  Model  with  Block  Sampling 
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phut  of  the  correlation  function  between  x(t)  end  y(t).  The 
relative  phase  0  between  the  input  and  reference  sinusoids  is 
chosen  to  be  zero  without  loss  of  generality  since  any  value  of  0 
can  be  represented  by  an  equivalent  value  of  time  delay  t. 

If  the  pre-envelopes  of  the  clipped  sinusoidal  signals  are 
represented  as  in  Chapter  III  by  the  Fourier  Series  rep re sent at i on r 


».  r2,  /  ,  J(-l)n  (2n  +  1)  wrt 

»x(t [  T^irfiT  * 

n=0 


(5.o; 


and 


y(t)  ■  7  I 

q=0 


IT  e 


-J(-Dq  (2q  +  1)  wRt 


(5-9) 


then,  from  Figure  5.8,  the  correlation  function  Zg(0  is 


M  N 

z6(t)  '  r  T  h  T  Px(mTP  +  \)  y(mTP  +  s +  t) 

m-l  i=l 


M  N  w  oo 

a.  i.  y  V  v  v  M)n+q 

2  MN  /  /  /  /  (2n  +  l)(2q  +  1) 

n  m=l  i=l  n=0  q=0 


j(-l)n  (2n  +  1)  UjOnTp  +  t±)  -j(-Dq  (2q  +  l)  ^(mT+t^-r) 


(5.10) 


where  n  denotes  the  input  odd  harmonics  and  q  denotes  the  reference 
odd  harmonics.  Changing  the  order  of  summation  of  Equation  (5.10)  gives 
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“  »  n+n  -j(-l)q  (2q+l) 

z6<t>  ■  %  y  y  (ri,Wi)(2q+i) — 

11  n=0  q=0 

1 
N 

1 
M 

For  specific  values  of  f_  and  fT,  Equation  (5. 11)  gives  the 
deterministic  correlation  function  as  a  function  of  time  delay  t. 

This  equation  shows  that,  in  general,  zg(T)  contains  contributions 
from  all  input  and  reference  odd  harmonics.  Each  of  these  terms  has 
an  amplitude  and  a  phase  which  are  functions  of  the  harmonic  order 
(n  and  q),  the  specific  frequencies  (fR  and  fj),  the  nonuniform 
sampling  instants  (t^),  the  block  sampling  period  (T^),  the  number 
of  samples  per  block  (N),  and  the  number  of  blocks  in  a  signal  length 
(M).  The  real  part  of  Zg(x),  ^(t),  which  represents  the  "in-phase" 
correlation  function,  thus  consists  of  the  sum  of  an  infinite  number  of 
sinusoidal  terms  whose  amplitudes,  frequencies,  and  phases  are 
determined  by  these  parameters.  As  will  be  shown  later  by  a  specific 
example,  this  results  in  a  correlation  function  with  a  "noisy"  carrier; 
i.e.,  the  basic  sinusoidal  carrier  term  at  frequency  fR(k  *  n  =  0)  in 
Equation  (5. 11)  is  modified  by  the  other  odd  harmonic  terms  to  create  a 
resultant  correlation  function  which  is  deterministic  but  irregular 


-J2itt.[(-l)q(2q+l)  fR-(-l)n(2n+l)  f  ] 
L=1 


&  -J27rmT  [(-l)q(2q+l)  f  - (-l)n(2n+l)  f  ].  (5. 11) 

1  6 
a=l 


in  shape. 
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Although  Equation  (5. 11)  is  exact  for  all  frequencies  f^  and 
fp,  it  is  quite  cumbersome  to  work  with,  in  general,  because  of  the 
four  summations.  However,  a  considerable  simplification  results  for 
the  important  Worst  Case  frequencies  described  by  Equations  (5-M> 
(5.5),  and  (5.6);  i.e.,  when 


and 


(*\ 


kfi  *  "iV 

tfR  =  mR  fp 


fR>  “  mdV 


(5.12) 

where  m^,  m^,  and  m^  are  integers.  This  simplification  results 

because,  when  these  three  conditions  are  satisfied,  the  term 
[(-l)q(2q  +  1)  fR  -  (-1)°  (2n  +  1)  fj],  which  represents  the  frequency 
difference  between  the  (2q  +  l)  reference  harmonic  and  the  (2n  +  l) 
input  harmonic,  reduces  to  a  multiple  of  i.e.,  whenever  the 

conditions  of  Equation  (5-12)  are  satisfied,  it  can  be  shown  that 


[(-l)q  (2q  +  1)  fR  -  (-1)"  (2n  +  1)  fj]  =  k  f  = 

P  P 

(5-13) 

where  k  is  a  positive  or  negative  integer  which  depends  on  ,  mp, 
and  n^.  When  the  Worst  Case  conditions  described  by  Equations  (5.12) 
and  (5.13)  exist,  the  final  summation  over  m  in  Equation  (5.11) 


reduces  to 
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!  V-  -J2«mTp[(-l)q  (2q+l)  fR  -  (-l)n  (2n+l)  fj] 

M  /_  e 
m=l 


M 

=  1  y  .  !  (MU) 

m=l 


for  all  k  and  n,  and  the  summation  over  i  can  be  expressed  in 
terms  of  the  sampling  coefficients  of  Chapter  IV  as 


x  'T’  -d2«ti[(-l)q  (2q+l)  fR  -  (-l)n  (2n+l)  fj] 

!i  L  * 

i=l 

N  -J2 k 

■  }  i  •  p  “  v 

i=l 


(5.15) 


nip  m 

where  k  is  an  integer  given  by  k  =  [(-l)q(2q+l)jp  -  (-l)n(2n+l)j~] . 

Thus,  for  the  Worst  Case  frequencies,  described  by  Equation  (5-12), 
Equation  (5. 11)  reduces  to 


B_  f  f  \  -J(-l)'Wl)  V. 

Z6^T'  "2  ./  L  (2n+l)(2q+l)  * 

n=0  q=0 


(5.16) 


Equation  (5.16)  shows  that  the  correlation  function  for  this  case  can 
be  written  in  terms  of  the  odd  input  and  inference  harmonics  and  of 
the  nonuniform  sampling  coefficients  a^,  discussed  in  Chapter  IV, 
thus  illustrating  the  important  role  of  these  coefficients.  The 
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and  whose  phases  are  determined  by  the  phase  of  the  complex  coefficient 


ak(q  n)'  ®^nce»  Chapter  IV,  the  a^  terms  can  vary  greatly  in  both 
amplitude  and  phase  from  term  to  term,  the  summation  given  by 
Equation  (5. IT),  in  general,  leads  to  an  irregularly  shaped  waveform 
which  is  periodic  at  the  period  of  the  reference  frequency  fR;  i.e., 
which  has  a  period  of  l/fR . 

Because  of  the  amplitude  weighting  of  unity,  the  most 
significant  term  in  Equation  (5 .17)  is  the  a^Q  ^  term  at  the 
fundamental  frequency  fD.  One  point  of  view  might  be  to  consider 
the  a^0  e  term  of  Equation  (5.17)  as  the  principal  carrier 

term  and  the  remaining  terms  as  additive  "noise"  terms,  although  for 
some  small  values  of  a^Q  the  noise  terms  may  be  larger  than  the 
principal  term.  Adopting  this  point  of  view,  if  we  wish  to  minimize 
the  odd  harmonic  correlation  responses,  we  should  try  to  keep  \(q  0) 
uniformly  low  in  some  sense  for  all  the  Worst  Case  frequency  conditions. 
The  a^Q  coefficient  is  a  single  efficient  in  the  block  sampling 


spectrum  where  the  integer  k 


(0,0) 


is  given  by 


k(0,0) 


^fR  " 


(5.18) 


If  a  number  of  reference  frequencies  an  input  frequencies  are  to  be 
considered,  Equation  (5.18)  can  be  used  to  compute  all  applicable 
values  of  k. 

One  approach  to  finding  an  optimum  set  of  a^Q  ^  coefficients 
would  be  to  try  to  choose  the  sampling  times  to  minimize  the 
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variance  of 


over  the  k's  given  by  Equation  (5. 18).  In 

,2 


over 


Chapter  IV,  a  technique  for  minimizing  the  variance  of  J  a^ 
an  entire  period  of  the  sampling  spectrum  was  developed.  The  problem 

I  1  2 

of  minimizing  the  variance  of  a^  |  over  an  arbitrary  set  of  k's 
is  a  considerably  more  complicated  problem  and  is  not  considered  in 
this  thesis.  Although  the  PRNG  sampling  technique  described  in 

I  I2 

Section  k.j  does  not  minimize  the  variance  of  |  av  |  in  a  localized 

2 

I  I 

k  region,  it  has  a  similar  property  in  that  the  values  of 


M 


are  suppressed  well  below  the  average  spectral  level  of  the  localized 
region  1  <  k  <  N/2.  If  the  differences  between  the  values  of  f 
and  fj  in  Equation  (5.18)  are  small  so  that  all  k(0,0)  terms  of 
interest  occur  in  the  region  1  <  k(0,0)  <  N/2,  the  PRNG  sampling 
technique  is  effective  in  substantially  reducing  the  contribution  of 


-JV 


in  the  correlation  function 


the  principal  term,  a^^  q)  1 
expansion  of  Equation  (5.I7). 

At  the  same  time,  the  PRNG  sampling  spectrum  is  reasonably 
uniform  in  amplitude  for  large  k  so  that  no  higher  order  odd  harmonic 
terms  of  Equation  (5. 17)  will  dominate  and  create  an  unusually  large 
correlation  peak. 


5  .U  Statistical  Analysis 

Although  the  correlation  function  Zg(r)  of  Figure  5-8  for 
fixed  values  of  fR  and  fj  is  a  deterministic  function  given  by 
Equation  (5. 11),  it  is  sometimes  useful  to  consider  the  statistical 
problem  of  finding  the  average  properties  of  the  correlation  functions 
for  all  reference  frequencies  which  respond  to  a  given  input  frequency 
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as  a  function  of  the  average  properties  of  a^  developed  in 
Chapter  IV.  For  example,  let  the  fixed  input  frequency  fj  satisfy 
the  Worst  Case  condition  Ufj  =  f  of  Equation  (5.12)  and 
consider  all  reference  frequencies  which  satisfy  Worst  Case  conditions 
*  *R  fp  “d  (fj.  -  fp)  =  fp.  The  corresponding  correlation 
functions  as  given  by  Equation  (5.16)  are 


*6<t)  =  %  f  f 


(-l)n+q  ak  -j(-l)q  (2q+l)  wrt 


it  —  — 

n=0  q-0 


Since  the  a^'s  vary  for  each  set  of  Worst  Case  frequencies  f^ 
and  fR,  this  analysis  will  assume  that  the  a^'s  are  random 
variables  with  a  mean  and  mean  square  given  by  Equations  (^.20)  and 
(U.17)  as 


where  II  is  the  number  of  samples  per  block  interval  T^,  and  K 
is  the  assumed  period  of  the  sampling  spectrum  a^.  Choosing  1  -  0 
as  an  arbitrary  value  of  carrier  phase,  the  average  value  of 


zAi)  at  x  =  0  is 
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where  the  average  is  taken  over  all  Worst  Case  input  and  reference 
frequencies.  The  average  squared  value  of  Zg(0)  taken  over  the  same 
frequencies  is 
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(5-19) 


where  the  notation  a^  ^  is  used  to  denote  the  dependence  of  a^ 
on  the  indices  n  and  q. 


Consider  the  average  a^^  a*k(p  r)  in  Equation  (5*19) 
over  all  input  and  reference  frequencies  satisfying  Uf^.  -  m^.  f^ 

and  (fj  -  fR)  -  fp-  We  can  consider  all  possible  Worst  Case 
input  and  reference  frequencies  by  averaging  both  and  over 

one  sampling  spectrum  period  K;  i.e.,  by  considering 


\(q,n)  a  k(p,r ) 
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(5.20) 


Substituting  the  expression  for  a^  and  interchanging  the  order  of 
summation,  Equation  (5.20)  becomes 
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where  ^  =  £[(-l)q  (2q+l)  -  (-l)n  (2n+l)]  and 

I2  =  "  (-1-)r  (2rfl)]  are  both  integers  for  all 

values  of  q,  n,  p  and  r.  But, 
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and 


K  -j|S  md[(-l)q(2q+l)ti  -  (-l)P(2pM)  tfc] 
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1  whenever  [(-l^3  (2q+l)td  -  (-l)p(2p+l)  thJ  =  0,+l,+2,, 


0  otherwise. 


(5-23) 
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*k(q  n)  a*k(p  r)  of  Ecluation  (5*21-)  will  be  non- zero  only  for  those 
values  of  q,  n,  p,  r,  and  t^  which  simultaneously  make  both 
Equations  (5.22)  and  (5.23)  non-zero.  From  Equation  (5-22)  and  (5.23), 
the  conditions  on  q,  n,  p,  r,  t^  and  t^  which  give  non-zero 
contributions  are? 


£[(-l)q  (2qtl)  -  (-l)n  (2n+l ) ]  t±  -  £[(-l)P  (2ptl)  -  (-l)r  (2r+l)]  th 


=  0 ,  +-1 ,  +2 ,  ... 


(5.2M 


L ( -1 )Q  (2q+l))  t±  -  [(-1)P  (2p+l ) ]  th  =  0,  +1,  +2 . 


(5.25) 


The  principal  non-zero  a^q  ^  R*k(p  r)  tenna  occur  when 


Q  -  P. 


n  =  r 


ti  '  V 


(5-26) 


thus  satisfying  both  Equations  (5.2l*)  and  (5.25)  by  making  them  equal 
to  zero.  For  these  terms, 

_  N  N 

Vq.n)  a*k(p,r)  -3  V  y  1  -  I  .  'for  ,  =  p.  h  =  r) 

i=l  h=l 


(5-27) 
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from  Equation  (5.21).  In  addition  to  the  principal  non-zero  terms 

given  by  the  Equations  (5.26)  and  (5.27)  which  always  exist  for  any 

set  of  sampling  times,  for  a  specific  set  of  £t^}  ,  a  few  other 

isolated  values  of  q,  p,  n,  r,  t^  and  th  may  satisfy  both 

Equations  (5-2l)  and  (5.25)  and  produce  small  but  non-zero  values  of 

V.  a*.  ,  v.  For  most  sets  of  sampling  times  <  t.  )  with  a 
q,n)  k(p,r)  *  \  U 

reasonable  degree  of  randomness,  the  effects  of  these  isolated 
contributions  are  second  order  and  will  thus  be  neglected  in  this 
analysis. 

Considering  only  the  principal  terms  given  by  Equation  (5.27), 
Equation  (5.I9)  reduces  to 
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(5-28) 


Equations  (5-18)  and  (5.28)  show  that  for  a  fixed  (arbitrary) 
value  of  x,  the  average  value  01  the  Worst  Case  correlation  function 
is  zero  and  their  mean  squared  value  is  equal  to  ^  independent  of 
specific  sampling  times  t^  if  the  averages  are  taken  over  a  complete 

sampling  spectrum  period  K.  Since  Zg(i)  is  a  complex  function 
composed  of  the  in-phase  and  quadrature  correlation  functions  $g(t) 
and  $g(t)  (see  Equations  2.1  -  2.6),  zg(0)  -  0  implies  than  both 

the  in-phase  and  quadrature  functions  also  have  a  zero  mean.  The  mean 
squared  value  of  l/N  given  by  Equation  (5.28)  is  the  average  value 
of  the  squared  "envelope1'  of  Zg(0)  since  by  Equation  (2.6)  |  ( T ) | 

is  defined  as  the  envelope  of  Zg(t)  in  the  usual  sense  of  the  rms 


158 


of  the  real  and  imaginary  parts  of  zg(i).  Equations  (5*18)  and  (5*28) 

thus  give  the  correlation  function  statistics  for  a  fixed  value  of  t 

for  all  possible  Worst  Case  frequency  combinations.  If  only  a  limited 

number  of  reference  frequencies  are  considered,  the  value  of 
2 

Zg(0)  -  l/N  given  by  Equation  (5.28)  may  not  be  representative 

of  these  reference  frequencies.  For  example,  if  PRNG  sampling  is  used 

and  if  the  principal  a,  terms  of  Equation  (5-19)  occur  in  the  small 

*  2 
k  suppression  region  shown  in  Figure  4.13,  the  value  of  Zg(o) 

can  be  considerably  less  than  l/N  as  will  be  shown  in  Section  5*5. 

Since  the  function  Zg(x)  is  composed  of  the  sum  of  a  large 

number  of  contributions  as  shown  in  Equations  (5*11)  and  (5.16),  one 

would  suspect  intuitively  that  its  measured  probability  distribution 

would  be  approximately  Gaussian.  This  was  experimentally  shown  to  be 

true  by  programming  a  digital  computer  to  simulate  the  correlation 

function  model  of  Figure  5-8  and  then  by  measuring  the  correlation 

function  statistics  for  the  Worst  Case  input  and  reference  frequencies. 

The  results  of  these  measurements  are  summarized  in  Section  5-5  for  an 

example  using  PRNG  sampling  with  N  =  98  samples  per  block  interval. 

Although  <t>g(0)  was  experimentally  shown  to  be  approximately 

Gaussian,  the  correlation  function  considered  as  a  function  of  t 

does  not  meet  the  requirements  of  narrow-band  Gaussian  noise  [13] 

because  its  spectrum  is  not  confined  to  a  narrow  band  due  to  the 

presence  of  the  odd  harmonic  response  terms.  As  a  result,  the 

correlation  function  $g(i)  cannot  be  approximated  by  the  usual 

narrow-band  assumptions  as  a  sinusoidal  function  with  a  slowly  varying 
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envelope  as  is  evident  from  an  examination  of  Equation  (5.16)  and  of 
the  example  of  Figure  5.9  in  Section  5.5.  Davenport  and  Root  [13] 
point  out,  however,  that  the  assumption  of  a  narrow-band  spectrum  is 
not  necessary  to  show  that  if  <tg(0)  and  $g(0)  have  Gaussian 
distri  ;ions,  the  function  defined  by 


6.(0)  ; 


zAo) 


-  V  (»6(0)} 


(»6(o)}2 


has  a  Rayleigh  probability  distribution 
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p(£)  =  ?  e 
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£.2  ° 


whose  mean,  mean-square,  and  variance  are  given  by 


t-  § 
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and 


2  a / 


<2-f)  ' 


where  og  is  the  variance  of  the  Gaussian  random  processes  <*g(0) 
or  5^(0).  Thus,  if  one  samples  the  outputs  ®g(t)  of  all  correlation 
functions  which  satisfy  the  Worst  Case  conditions  of  Equation  (5 .12) 
at  any  fixed  time  t,  those  samples  will  have  an  approximate  Gaussian 
distribution  with  a  zero  mean  and  a  variance  of  ^  .  Likewise, 
samples  of  the  output  "envelopes"  |zg(-r)|  ^or  the  same  channels  taken 
at  any  fixed  time  t  will  have  an  approximate  Rayleigh  distribution 


with  a  mean  of  Vfe  and  a  variance  of  [(2  -  — )  •  This 

approximate  statistical  model  of  the  Worst  Case  correlation  function 
samples  was  verified  experimentally  as  will  be  shown  in  Section  5-5* 
Although  the  statistical  model  of  the  correlation  function  as 
a  wide-band  Gaussian  process  with  a  Rayleigh  envelope  distribution  is 
seen  to  fit  the  system  model  of  Figure  5.1,  the  concept  of  zg(t)  as 
an  envelope  normally  has  real  significance  only  for  narrow-band 
processes.  For  example,  if  it  is  desired  to  compute  the  peak  value  of 
^(t)  over  one  period  (see  Figure  5*9)»  it  might  be  necessary  to 
redefine  the  envelope  of  in  a  more  meaningful  way. 

Equation  (5.28)  shows  that  if  we  relate  the  average  properties 
of  the  correlation  functions  for  all  Worst  Case  frequencies  to  the 
average  properties  of  the  sampling  spectrum  coefficients  from  Chapter  IV, 
the  mean-squared  value  of  Zg(0)  is  given  by 

2 

z6(0)  =  1/N  , 

where  N  is  the  number  of  samples  per  block.  We  now  wish  to  extend 
this  analysis  to  the  non-worst  case  frequency  combinations  shown  in 
Figure  5.6  and  5.7  and  discussed  qualitatively  in  Section  5.2  in  order 
to  obtain  a  first  order  comparison  between  the  correlation  function 
statistics  in  these  various  conditions.  For  signals  whose  duration  T 
is  long  relative  to  a  block  length  T^,  the  key  to  the  differences  in 
correlation  for  the  different  frequency  combinations  will  be  shown  to 
lie  in  the  simulation  over  m  in  Equation  (5. 11);  i.e.,  in 
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l6l 


A  -j2*mT  [(-l)q(2q+l)  fR  -  (-l)n(2n+l)  f  ] 

,:5  [  ' 

m=l 

where 

M  =  number  of  sampling  blocks  of  length  in  one  signal 
duration  T, 
fj  =  input  frequency 
and 

fR  =  reference  frequency. 

It  was  shown  in  Equation  (5.1*0  that  F  =  1  for  all  q  and 
n  for  the  Worst  Case  frequencies  because  the  product 

P  =  T  [(-l)q(2q+l)  f.,  -  (-l)n(2n+l)  fT]  is  always  equal  to  an 
p  K  1 

integer  for  this  case.  If  either  fR  or  f^  do  not  satisfy  the 
Worst  Case  conditions ,  P  will  not  be  an  integer  for  all  values  of 
k  and  n  and  the  summation  F  will  be  less  than  unity  in  absolute 
value,  thus  acting  as  an  attenuation  factor  in  the  correlation  function 
of  Equation  (5- 11).  For  large  M,  the  attenuation  factor  F,  whose 
average  value  is  proportional  to  l/M,  will  be  a  significant  smoothing 
factor. 

As  an  example  of  the  effect  of  the  attenuation  factor  F, 
consider  the  third  harmonic  response  frequencies  shown  in  Figure  5.6; 
i.e.,  those  satisfying 
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and 
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where  i  =  1  or  2 .  From  Equation  (5.29),  F  becomes: 


(5-30) 
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=  M  /  * 
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and 


F 


=  1  when  (2q+l)  *  3,  9,  15,  21,  etc, 


<  pj  otherwise  if  M  >  3. 


If  we  assume  M  is  large  and  neglect  all  values  of  n  and  q  for 
which  F  |  <  j-j  ,  the  correlation  function  at  t  =  0  for  the 
frequencies  satisfying  the  conditions  of  Equation  (5.30)  can  be 
written  [from  Equation  (5. 11)]  as: 

8  r1  V1  (-l)n+<5  *».(-  n\ 

z6(°)  =  *2  I  1  (2n*l)W?l  • 

n*0  q=0 

and  the  mean  squared  value  becomes 
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which  is  1/9  of  the  Worst  Case  value  given  by  Equation  (5.28). 
Under  the  sane  assumptions,  it  can  be  shown  that  for  the  fifth 
harmonic  response  frequencies  of  Figure  5*6 

I  , ,  2  1  1 

h(o)  *  55  ■»  ■ 


and,  in  general,  for  the 


Kth 


harmonic  response 


z6{0) 


1 

N  * 


Thus,  for  the  situation  represented  by  Figure  5.6  in  which  the  input 
frequency  fj  satisfies  4fj  -  m^f^,  the  output  power  for  the 
correlation  function  corresponding  to  the  Kth  harmonic  is  inversely 
proportional  to  K2  as  a  first  approximation  for  large  M. 

Next,  consider  the  conditions  represented  by  Figure  5*7  where 
fj  does  not  satisfy  the  Worst  Case  input  frequency  condition;  i.e., 
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where  4f^  /  m^f  .  ^rie  referenre  frequencies  with  the  maximum 

response  are  the  ones  for  which  fp  =  f  j  +  .  For  these 

frequencies,  F  from  Equation  (5-29)  is 


x  Y  -j2JtmT^[(.l)q(2q+l)fI-(.l)n(2n+l)fI+ (-l)q  (2q+l)  mdfp] 
F  =  M  /  e 


i  -j2irmT  [(-l)q(2q+l)  f  -  (-l)n(2n+l)  f  ] 

m  L  6 

m-1 


-  1  when  q  =  n 


<  1  otherwise. 


Although  j F  J  can  occasionally  be  nearly  equal  to  unity  for  q  /  n, 

we  can  neglect  these  terms  with  the  result  that  our  estimate  of 
2 

zg(0)  will  be  noisy;  i.e.,  have  a  larger  variance.  Neglecting 
all  terms  for  which  F  /  1,  the  complex  correlation  function  at 
t  -  0  for  this  case  becomes: 


(0)  =~  8  t  \ln4l 

«  L,  (2n+l)2 


and  its  mean  squared  value  becomes: 
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(5-31) 


which  can  be  compared  with  a  value  of  l/N  for  the  Worst  Case 
condition  of  Figure  5*6.  Therefore,  if  the  input  frequency  does  not 
satisfy  the  condition  Uf^  =  m^f^,  samples  of  «g(0)  from  all 

reference  channels  satisfying  fD  -  fT  +  m.f  approximately  satisfy 

n  1  Q  p 

a  Gaussian  distribution  with  a  zero  mean  and  a  variance  of  , 

3N 

which  is  2/3  of  the  Worst  Case  variance  of  ^  as  shown  in 
Figure  5.6.  Likewise,  the  envelope  of  Zg(0)  defined  as  |z^(0 )| 
has  an  approximate  Rayleigh  distribution  with  a  mean  of 

(VT  *V¥)  a  variance  of  [(2.|).^]. 

Using  a  similar  analysis,  all  the  third  harmonic  response 

frequencies  shown  in  Figure  5.7  can  be  shown  to  have  a  mean  squared 

value  of  1/9  of  the  maximum  response  channel.  Likewise,  the 

variance  in  the  fifth  harmonic  channels  is  approximately  ~ 

and,  in  general,  for  the  Kth  harmonic  response  the  variance  is 

approximately  given  by  “2  ’  * 

K 

Although  the  statistical  analysis  in  this  section  is  not 


rigorous  in  the  sense  that  the  correlation  functions  it  represents 
are  actually  deterministic,  it  has  utility  in  that,  for  long  signals 
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with  M  »  1,  it  gives  an  approximate  but  useful  model  of  the  average 
correlation  and  its  variation  for  all  Worst  Case  frequencies.  For 
signals  whose  time  duration  T  is  not  long  relative  to  a  block 
interval  Tp,  the  effects  of  the  finite  T  often  dominate  the  higher 
order  harmonic  effects  as  will  be  shown  by  some  examples  in  Section  5-5- 

5-5  Example  Using  Pseudorandom  Block  Sampling 

In  Section  U.7,  a  method  for  generating  a  repetitive  set  of 
non-uniform  sampling  times  using  pseudorandom  number  generators 
(FRNG's)  was  described.  Although  this  sampling  scheme  was  shown  to  be 
non  minimum  variance  in  general,  it  does  possess  a  high  degree  of 
randomness ,  has  a  unique  spectral  amplitude  suppression  characteristic 
for  small  k  and  has  a  good  deal  of  flexibility  in  the  selection  of 
the  number  of  sampling  times  per  block  N,  and  the  number  of  possible 
sampling  instants  per  block  K.  A  digital  computer  simulation  of  a 
quadrature  correlation  function  for  clipped  sinusoids  with  pseudo¬ 
random  block  samp!  ing  was  made  in  order  to  experimentally  verify  the 
theoretical  model  of  Figure  5.8.  A  summary  of  the  results  of  that 
simulation  are  presented  in  this  section. 

The  block  sampling  parameters  chosen  for  the  simulation  are 
summarized  in  the  following  table: 
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Simulation 

Value 

Parameter 

Example  1 

Example  2 

Sample  Pairs  Per  Block,  N 

96 

2k 

Possible  Sampling  Instants  Per  Block,  K 

12,  192 

ZOkQ 

Blocks  Per  Signal  Duration,  M 

3 

3 

PRNG  Codes 

Several 

Several 

As  previously  mentioned,  a  typical  correlation  function  for 
clipped  sinusoids  and  pseudo-random  sampling  has  a  noiselike 
characteristic  because  of  the  presence  of  odd  harmonic  components  and 
tends  to  lose  the  narrowband  appearance  of  a  linear  correlation 
function  with  uniform  sampling.  This  fact  is  illustrated  in  Figure  5.9 
where  one  period  of  a  typical  correlation  function  from  the  simulation 
program  for  Example  1  is  shown.  The  two  quadrature  phases  of  the 
correlation  function,  shown  by  the  solid  curves  are  seen  to  be  n/2 
phase  shifted  translations  of  each  other.  Although  the  fundamental 
2«  "carrier"  periodicity  is  still  evident  in  these  functions,  their 
"noisy"  shapes  vary  considerably  from  the  conventional  sinusoidal 
carriers  found  in  linear  correlation  functions.  The  dashed  curve  in 
Figure  5*9  shows  that  the  narrow  band  definition  of  envelope,  as  the 
square  root  of  the  sum  of  the  squares  of  the  two  quadrature  components, 
does  not  have  the  usual  geometrical  interpretation  here  of  a  slowly 
varying  function  passing  through  the  peaks  of  the  carrier. 

Some  examples  of  correlator  output  functions  are  shown  in 
Figure  5-10  for  Example  2  (N  =  24),  where  't'g(i)  is  plotted  for 
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one  output  repetition  period.  The  typical  Worst  Case  correlation 
function  of  Part  (a)  shows  a  high  degree  of  randomness  in  the  waveform 
similar  to  that  of  Figure  5.9.  The  correlation  function  shown  in 
Part  (b)  is  for  a  perfectly  matched  frequency  case  in  which  the 
reference  frequency  fR  is  equal  to  the  input  frequency  f ^ .  For  a 
correlation  function  with  uniform  sampling,  the  plot  of  Part  (b) 
would  have  a  linear  sawtooth  characteristic  of  period  2n. 

The  deviation  lrom  linearity  (i.e.,  the  granularity)  of  the  correlation 
function  of  Figure  5.10(b)  as  well  as  the  stepwise  character  of 
Figure  5.10(a)  is  due  both  to  the  nonuniformity  of  the  sampling  function 
and  to  the  relatively  small  number  of  samples  per  block  N. 

The  computer  simulation  consisted  of  computing  the  exact 
correlation  function  by  simulating  two  square  waves  with  frequencies 
fR  and  fj,  time  shifting  the  reference  square  wave,  sampling  both 
signals  at  the  pseudo-random  sampling  times,  and  correlating  the 
sampled  signals.  The  curves  of  Figures  5-9  and  5-10  were  generated 
in  this  way  using  100  increments  of  urt  in  a  2n  output  interval. 

In  addition,  a  second  computer  program,  useful  for  studying  the  effects 
of  specific  harmonics,  was  written  to  compute  the  correlation  function 
term  by  term  from  Equation  (5.16).  Figure  5-ii>  which  shows  a  typical 
Worst  Case  output  correlation  function  from  the  small  k  suppression 
region  for  Example  1  (N  -  96),  shows  a  four-step  comparison  of  the 
simulation  model  with  the  Equation  (5.16)  expansion.  The  dashed 
curve,  which  is  identical  in  all  three  plots,  represents  the  exact 
output  correlation  function  while  the  four  solid  curves  represent 
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approximations  using  Equation  (5.1 6)  with  1,  5>  10,  and  55  terms. 

The  0)  term  shown  in  the  top  plot  is  not  a  major  contributor 

since  its  peak  amplitude  is  only  about  1.6  percent.  This  is 
characteristic  of  all  Worst  Case  correlation  functions  in  the  small 
k  region.  The  five-term  approximation  which  contains  the  fundamental, 
third,  and  fifth  harmonic  terms  is  seen  to  follow  the  general  shape  of 
the  dashed  curve,  but  the  approximate  correlation  is  inaccurate  by  as 
much  as  6  percent  and  the  high  frequency  variations  are  not  present. 

The  10-term  approximation,  which  has  reduced  the  maximum  error  to 
about  4  percent,  is  starting  to  conform  a  little  more  closely  to  the 
exact  wave  shape.  The  55-term  approximation  shows  the  approximation 
slowly  converging  to  the  true  correlation  function.  This  relatively 
slow  convergence  shows  that  a  significant  number  of  the  a^  terms  in 
Equation  (5.16)  are  important  and  emphasizes  the  importance  of 


maintaining  a  low  variance  of 


*k 


over  the  entire  spectral  period 


while  keeping  the  a^g  q)  term  low  by  small  k  suppression  as 
discussed  in  Sections  4.6  and  4.7. 

The  significance  of  the  small  k  spectral  suppression  effect, 
a  characteristic  of  all  PRNG  generated  sequences,  is  seen  by  an 
examination  of  Figure  5-12  which  shows  three  plots  of  peak  output 
correlation  (over  a  =  2n  interval)  versus  reference  frequency 
for  fixed  input  frequency.  Figure  5.12(a)  shows  a  typical  example 
for  reference  frequencies  with  a. 


terms  in  the  high  k  region 

satisfj 

reference  frequencies  marked  are  the  computed  Worst  Case 


and  a  Worst  Case  input  frequency  satisfying  4f  ^  -  m^f^ ,  The 


Figure  5.12  Typical  Experimental  Plots  of  Peak  Correlation  Vs 
Reference  Frequency  for  Pseudo-Random  Block 
Sampling  (N  =  96,  M  =  3) 


17*» 


frequencies  satisfying  l*fR  <■  ^f  ,  Because  of  the  finite  signal 
duration  of  T  seconds,  the  typically  large  peak  correlations 
associated  with  the  Worst  Case  reference  frequencies  extend  over 
approximately  a  reference  frequency  region  of  width  l/T.  Although 
the  maximum  of  the  Ja<sged  peak  correlation  versus  reference  frequency 
curve  occurs  near  the  ^W^  points,  they  do  not  always  fall  exactly 
on  the  Worst  Case  frequencies  (e.g.,  see  (w^)  ).  Although  the  © 

frequencies  tend  to  have  the  largest  correlations,  sometimes,  as  Jhown 
at  ^W^  ,  the  odd  harmonic  terms  add  in  such  a  way  that  the  peak 
correlation  never  achieves  a  significant  value.  As  shown  in 
Figure  5 -12 (a),  the  higher  order  odd  harmonic  effects  occurring  between 
the  ®  frequencies  are  essentially  dominated  by  the  l/T  spreading 
effects  for  the  particular  parameters  (high  N,  low  M)  of  this 
example.  The  example  shown  in  Figure  5.12(b),  which  covers  reference 
frequencies  in  the  small  k  suppression  region,  shows  typical 
reductions  in  the  Worst  Case  frequency  peak  correlation  of  about 
33  percent  relative  to  those  for  Figure  5.12(a).  This  is  attributable 
to  the  absence  of  significant  0)  t,erms  in  this  r«gion»  The 

average  correlation  level  in  the  odd  harmonic  response  region  between 
Worst  Case  references,  as  expected,  did  not  change  from  Figure  5.12(a) 
to  Figure  5- 12(b)  because  the  terms  for  these  harmonics  are 

generally  found  in  the  high  k  spectral  region  in  both  cases.  The 
example  of  Figure  5 .12(c)  for  a  non-Worst  Case  input  frequency  shows 
much  less  fluctuation  of  peak  correlation  versus  reference  frequency 
because  the  Wirst  Case  frequencies  are  substantially  reduced  and  the 
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odd  harmonics  are  more  randomly  spaced  in  the  reference  frequency. 

The  odd  harmonic  effects  predicted  in  Section  5*^  would  be  more 
prominent  for  examples  with  longer  T  and  larger  M. 

Section  5.h  discussed  the  possibility  of  assuming  the  i  0 
samples  of  all  Worst  Case  correlation  functions  >t>g(0)  have  a  Gaussian 
probability  distribution  with  a  variance  of  l/2N.  This  assumption 
was  verified  experimentally  as  shown  in  Figure  5*1-3  using  the  computer 
simulation  program.  Figure  5* 13(a)  shows  the  results  of  300 
measurements  for  Example  1  (N  -  96)  taken  over  a  wide  range  of  Worst 
Case  reference  frequencies,  all  with  a^Q  ^  terms  in  the  high  k 
region.  The  experimental  histogram  shown  by  the  solid  curve  had  a 
measured  standard  deviation  of  0  -  7.5  percent  or  equivalently  a 

2  1  Oft 

variance  of  a  -  O.OO5625  =  which  deviates  from  the 

theoretical  value  of  ~  by  only  8  percent.  The  histogram  is  seen  to 
match  the  dashed  Gaussian  distribution  with  a  theoretical  standard 
deviation  of  7*21  percent  very  closely  except  for  very  low  correlations 
around  zero,  a  problem  which  is  believed  to  be  due  to  the  fact  that 
the  granularity  of  the  correlation  function  is  roughly  equal  to  the 
cell  size  used  for  plotting  the  histogram.  The  absolute  value  of  the 
correlation  was  used  instead  of  its  full  positive  and  negative  range 
in  order  to  double  the  sample  size  in  each  cell.  Figure  5  13(b)  shows 
a  plot  of  the  histogram  of  the  "envelope''  function  Zg(0)  for  the 
same  reference  frequencies  used  in  Figure  5.13(a)  with  a  sample  size  of 
150.  The  measured  square  value  was  Zg(0) 


probability  density  probability  density 
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Figure  5-12  Measured  Correlation  Histograms  for  Worst  Case 

Frequencies  in  Large  k  Region  -  FRNG  Sampling  (N  -  96) 


compared  with  a  theoretical  value  of  l/N  =  0.0104.  The  data  is 
seen  to  fit  fairly  well  to  the  theoretical  Rayleigh  distribution  with 
£2  =  l/N  =  0.0104. 

Figure  5*44  contains  experimentally  measured  histograms  for 

Worst  Case  reference  frequencies  from  the  small  k  spectral  suppression 

region.  Since,  for  the  high  k  region,  Equation  (5-28)  shows  that  the 

64  /• 

a^0  terms  contributed  about  -£■  -  65.5  percent  of  the  total 
variance  measurement,  Figure  5-44  shows  that  the  a^Q  q)  contribution 
has  been  reduced  by  about  90  percent  in  the  low  k  region. 

Because  of  the  noisy  nature  of  the  envelope  function  zg(^)|  > 
the  addition  statistic  MAX  ,  the  maximum  value  of  the  output 

correlation  function  over  cr  "it  period  was  measured  and  plotted  for 
571  Worst  Case  channels  in  the  small  k  region  as  shown  in  Figure  5.I5 
for  samples  of  data  from  Example  1  (N  =  9^)*  A  comparison  of  +he  peak 
statistic  histogram  of  Figure  5.I5  with  the  histogram  of  |zg(0)  in 
Figure  5.14  shows  that  the  mean  of  the  peak  distribution  is  about 
twice  as  large  as  the  mean  of  the  |zg(0)|  distribution. 
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(b) 

Figure  5-i^  Measured  Correlation  Histograms  for  Worst  Case 
references  in  Small  k  Suppression  Region  - 
PRNG  Sampling  (N  -  96) 
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Figure  5.15  Correlation  Peak  Histograms  for  Worst  Case 
References  in  Small  k  Suppression  Region 


CHAPTER  VI 


SUMMARY  AND  CONCLUSIONS 

6.1  Statement  of  the  Problem 

The  general  topic  of  the  research  reported  in  this  thesis  is  the 
study  of  nonuniform  block  sampling  theory  and  its  application  to  the 
study  of  correlation  functions  of  sampled  signals  that  are  limited  to 
a  frequency  band.  The  research  evolved  from  the  specific  problem  of 
trying  to  reduce  the  odd  harmonic  terms  which  appear  in  the  correlation 
function  of  hard  limited  and  uniformly  sampled  sinusoidal  signals. 

These  odd  harmonic  terms  arise  because  of  the  synchronism  between  the 
periodicity  of  the  uniform  samples  and  the  periodicity  of  the  odd 
signal  harmonics  produced  by  hard  limiting. 

6.2  Approach 

The  general  approach  to  the  problem  was  the  use  of  nonuniform 
sampling  to  reduce  the  undesirable  periodicity  effects  and  thus  to 
reduce  the  correlation  function  odd  harmonic  terms.  A  particularly 
useful  class  of  nonuniform  sampling  functions  for  this  application  was 
found  to  be  block  sampling  functions;’1  i.e  ,  short  sequences  (or 
blocks)  of  nonuniform  samples  which  are  repeated  periodically. 

As  a  preliminary  to  a  study  of  block  sampling  functions  and 
their  application  to  the  odd  harmonic  problem,  a  detailed  study  of  the 


general  properties  of  correlation  functions  for  large  TW  signals 
with  uniform  sampling  was  made.  These  results  were  then  extended  to 
the  case  of  sinusoidal  signals  with  uniform  sampling  to  illustrate  the 
role  of  both  clipping  and  sampling  in  the  generation  of  the  undesirable 
odd  harmonic  terms  in  the  correlation  function.  The  uniform  sampling 
correlation  function  results  along  with  the  results  of  the  block 
sampling  spectrum  analysis  then  formed  the  basis  for  analysis  of 
correlation  functions  of  sinusoidal  signals  with  nonuniform  block 
sampling. 

6.3  Discussion  of  Results 

In  the  study  of  correlation  functions,  two  alternative 
definitions  of  sampled  correlation  functions  were  compared,  the 
regular  correlation  function  given  by: 

M 

z(mts)  =  i  J  *<»»,)  y("T.) 

m-1 

and  the  qiadrature  correlation  function  given  by: 

M/2 

zq(MTs)  a  k  )  W2mTg)  y(2mTj  4  ^(2mTg)  y(2mT.)], 
m=l 

where  x(t)  and  y(t)  are  the  signals  of  duration  T  to  be  correlated 
and  x(t)  and  y(t)  are  oheir  Hilbert  transforms.  The  regular 
correlation  function  requires  M  uniform  signal  samples  taken  at  T 

s 

second  intervals  in  T  seconds  to  compute  each  point  of  z(MT  ), 

s 


whereas  the  quadrature  correlation  function  requires  M/2  sample 
pairs  of  the  signal  and  its  Hilbert  transform  taken  at  2T,  second 

o 

intervals  to  compute  z^(MTs).  It  was  shown  that  the  basic  difference 
between  the  two  correlation  function  definitions  is  that  quadrature 
correlation  eliminates  the  second  harmonic  in  the  spectrum  of  the 
product  x(t)  •  y(t).  One  implication  of  the  second  harmonic 
elimination  for  unclipped  signlas  was  that  (for  sampling  faster  than 
the  Shannon  rate)  the  quadrature  correlation  function  eliminates  second 
harmonic  sampling  errors  that  occur  m  the  regular  correlation  function. 
Another  implication  is  that  the  sampled  quadrature  correlation  function 
variance  for  narrowband  stationary  random  signals  is  independent  of 
sampling  rate,  whereas  the  regular  correlation  function  variance 
versus  sampling  rate  response  contains  large  peaks  at  some  values  of 
fg  due  to  the  presence  of  the  second  harmonic.  Thus,  for  no  increase 
in  the  number  of  computations  per  correlation  function  point,  quadrature 
correlation,  in  general,  achieves  more  accurate  correlation  for 
deterministic  signals  and  lower  variance  at  all  sampling  rates  for 
narrowband  random  process  signals. 

A  comparison  of  the  sampled  quadrature  correlation  function 
variance  for  clipped  and  unclipped  narrowband  random  processes  showed 
that,  although  some  second  harmonic  effects  appeared  for  the  clipped 
signal  case,  the  variance  could  be  reduced  by  as  much  as  3  dB  by 
sampling  faster  than  the  Shannon  rate. 

In  extending  the  study  of  quadrature  correlation  function 
properties  to  sinusoidal  signals,  the  effects  of  linear  sampling  and 
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clipping  were  derived  separately  and  the  combined  effects  of  both 
sampling  and  clipping  were  examined.  Uniform  sampling  was  shown  to 
make  the  correlation  function  periodic  (as  a  function  of  frequency) 
with  a  repetition  period  equal  to  the  sampling  rate,  while  clipping 
was  shown  to  introduce  extraneous  correlations  for  input  signals  with 
frequencies  which  are  odd  multiples  of  the  reference  signal  frequency. 
When  both  clipping  and  sampling  are  employed,  each  odd  harmonic 
correlation  response  is  repeated  at  periodic  intervals  which  are 
related  to  the  sampling  rate,  thus  giving  rise  to  the  odd  harmonic 
problem.  When  a  sinusoid  is  correlated  with  a  narrowband  random 
function,  the  quadrature  correlation  function  variance  versus  sampling 
rate  curve  approaches  its  high  sampling  rate  asymptote  more  slowly 
than  in  the  large  TW  signal  case,  and  the  second  harmonic  ripple  is 
more  pronounced. 

After  proposing  nonuniform  block  sampling  as  a  possible  solution 
to  the  odd  harmonic  problem  in  the  correlation  function  of  clipped  and 
sampled  sinusoids,  the  spectral  properties  of  block  sampling  were 
analyzed.  Although  other  studies  have  considered  block  sampling  from 
the  points  of  view  of  signal  reconstruction  from  sampled  values 
[4l,  2^],  sampling  interval  statistics  [7],  and  sampled-data  control 
system  analysis  [36],  they  did  not  analyze  the  basic  properties  of  the 
sampling  function  or  its  frequency  spectrum  which  are  necessary  for  a 
study  of  correlation  functions  of  nonuniformly  sampled  signals. 
Accordingly,  this  thesis  examined  in  detail  the  complex  frequency 
spectrum  properties  of  ideal  (zero  sampling  pulse  width)  block  sampling 


functions  including  amplitude  and  phase  characteristics  of  each  spectral 

component,  spectral  periodicity  conditions,  conservation  of  spectral 

energy  properties,  and  conditions  for  minimum  variance  of  the  sampling 

power  spectrum.  The  significant  sampling  spectrum  parameter  was  shown 

to  be  the  complex  sampling  coefficient  a^,  given  by 

t 

N  -J2jtksr 

a  -  i  f  e  » 

ak  N  /_  6  » 

n-1 


where  k  is  an  integer  designating  the  k  spectral  component,  t 

are  the  nonuniform  sampling  times,  is  the  block  repetition  period 

in  seconds,  and  N  is  the  number  of  samples  per  block.  Conditions 

on  tn/T  for  all  values  of  n  were  derived  such  that  a^  will  be 

periodic;  i.e.,  such  that  a^  =  a^4K  for  some  K  and  all  k.  A 

power  spectrum  conservation  property  was  established  stating  that  the 

2 

average  value  of  a^  taken  over  one  spectral  period  K  was 

equal  to  l/N  independent  of  the  specific  sampling  times  t  .  Thus, 
each  block  sampling  function  with  N  samples  per  block  (of  which 
uniform  sampling  is  a  special  case)  merely  redistributes  a  fixed  amount 


of  spectral  energy  as  a  function  of  frequency.  Although  the  average 
2 

value  of  a.  was  shown  to  be  a  constant  equal  to  l/N  independent 

2 

of  the  sampling  times  tn,  the  variance  of  a^  was  found  to  be 

very  dependent  on  the  specific  sampling  times  t  .  Accordingly,  a  set 

of  conditions  on  t  were  derived  which  would  minimize  the  variance  of 

n 


*k 


If  o. 


denotes  the  variance  of 


■k 


over  one  spectral 


period  K,  it  was  shown  that 


* 
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where  the  lower  limit  is  determined  by  the  minimum  variance  conditions 
and  the  upper  limit  is  set  by  uniform  sampling.  As  an  example  of 
sampling  spectrum  shaping;  i.e.,  of  controlling  the  values  of  a^ 
in  specific  spectral  regions,  by  choice  of  t  ,  a  block  sampling 
function  example  was  introduced  which  used  a  pseudo-random  number 
generator  (FRNG)  for  sampling  time  selection.  Although  this  sampling 
scheme  does  not  yield  a  minimum  variance  sampling  spectrum,  it  achieves 
a  unique  spectrum  shaping  by  suppressing  the  values  of  a^  for  small 
k,  a  property  that  is  important  in  the  reduction  of  odd  harmonic 
responses  in  correlation  functions  of  clipped  sinusoidal  signals. 

Finally,  having  derived  the  fundamental  properties  of  quadrature 
correlation  functions  and  of  nonuniform  block  sampling,  the  application 
of  block  sampling  to  the  reduction  of  odd  harmonics  in  correlation 
functions  of  clipped  sinusoids  was  studied.  Considering  the  quadrature 
correlation  function  z(t,  fj,  f^)  between  two  clipped  sinusoidal 
signals  with  block  sampling  as  a  function  of  the  two  sinusoidal 
frequencies  f^  and  f^,  it  was  shown  that,  although  the  odd  harmonic 
responses  are  greatly  reduced,  certain  combinations  of  input  frequency 
fT  and  reference  frequency  fD  result  in  correlations  whose  average 
value  is  greater  than  those  for  other  values  of  fj  and  f^.  These 
significant  combinations  of  and  fR,  which  are  called  Worst  Case 

combinations,  are  given  by  the  equations 
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where  ,  m^,  and  md  are  integers.  Based  on  an  analytical  model 
similar  to  one  used  by  Tou  [36],  an  expression  relating  the  correlation 
function  at  these  Worst  Case  frequencies  to  the  sampling  function 
coefficients  and  the  signal  odd  harmonics  was  obtained.  This 

expression  was  useful  for  studying  the  effect  of  each  odd  harmonic 
on  the  correlation  function  as  well  as  for  demonstrating  the 
effectiveness  of  the  small  k  suppression  property  of  FRNG  sampling 
functions  in  reducing  the  correlation  function  amplitude  for  the 
Worst  Case  conditions.  This  model  was  also  useful  for  obtaining  scr.ie 
average  correlation  characteristics  taken  over  all  Worst  Case  frequency 
combinations.  Some  experimental  results,  based  on  a  digital  computer 
simulation,  were  obtained  to  verify  the  theoretical  correlation  function 
model.  Both  the  analysis  and  the  experimental  results  confirmed  that, 
by  proper  choice  of  sampling  parameters,  block  sampling  could  reduce 
the  undesirable  odd  harmonic  correlation  function  responses  due  to 
clipping  and  sampling  by  any  desired  degree. 
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6 ,h  Suggestions  for  Further  Research 

The  quadrature  correlation  function  properties  obtained  in 
Chapters  II  and  III  were  derived  for  single  input  and  reference 
signals.  In  many  applications,  it  is  important  to  consider  sums  of 
signals  [23]  or  signals  plus  random  noise  processes  [30]  instead  of 
single  signals.  In  the  possible  extension  nf  the  results  of 
Chapters  II  and  III  to  these  cases,  of  particular  interest  would  be 
the  effects  of  hard  limiting  and  uniform  sampling  on  the  quadrature 
correlation  functions  of  narrowband  signals--especially  sinusoids. 

Throughout  this  thesis,  only  extreme  signal  quantization 
(i.e.,  hard  limiting)  was  considered.  However,  because  of  the  wide 
spread  use  of  digital  computers  and  other  digital  devices  to  calculate 
correlation  functions ,  it  is  also  important  to  study  other  levels  of 
quantization.  For  example,  since  all  computers  have  finite  word 
length,  some  degree  of  round-off  (i.e.,  quantization)  is  always 
present . 

In  the  analysis  of  the  sampling  function  spectrum  of  Chapter  IV, 

i  2 

the  sampling  spectrum  variance  o„  was  minimized  over  a  complete 
spectral  period.  Many  times,  however,  an  entire  spectral  period  need 
not  be  considered.  In  these  cases,  it  would  be  more  meaningful  to 
minimize  the  spectrum  variance  over  an  arbitrary  spectral  range  at  the 
expense  of  spectral  amplitudes  out  of  the  range  of  interest.  In 
addition  to  optimizing  the  spectrum  variance  by  choice  of  sampling 
times,  it  might  also  be  possible  to  fix  the  nonuniform  sampling  times 
and  minimize  the  sampling  spectrum  variance  by  the  choice  of  amplitude 


weighting  coefficients  for  each  sample.  Another  more  difficult  related 
problem  is  the  synthesis  problem  of  specifying  the  exact  desired  shape 
of  the  sampling  spectrum  and  finding  a  set  of  sample  times  which  yields 
this  spectrum  (or  a  closest  fit  to  it). 

In  Chapters  IV  and  V,  out  of  many  possible  specific  nonuniform 
block  sampling  schemes ,  only  PRNG  sampling  was  studied  in  detail . 

Many  other  intersting  sampling  schemes  for  further  study  suggest 
themselves,  such  as  using  a  probability  law  (Gaussian,  Poisson,  etc.) 
to  select  the  nonuniform  sampling  times  or  restricting  the  freedom  of 
each  sample  within  the  sampling  block.  For  the  PRNG  sampling  scheme 
studied,  additional  comparative  analysis  of  the  effects  of  specific 
pseudo- random  codes  are  suggested. 


’ 
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APPENDIX 


CORRELATION  FUNCTION  VARIANCE  FOR  RANDOM  PROCESS 
INPUT  SIGNALS  WITH  RECTANGULAR  POWER  SPECTRUM 


We  wish  to  compute  the  correlation  function  variance  for  a 
typical  example  where  the  input  and  reference  signals  are  assumed  to 
be  random  variables  with  Gaussian  distributions  and  to  be  represented 

by 

x(t)  =  an  n(t)  =  CTn(Ni  cos  “o1  "  N2  sin  (A.l) 

and 


y(t)  =  s ( t )  =  E^  cos  u)Qt  -  Eg  sin  wQt  .  (A. 2) 

Assume  that  the  power  spectrum  of  both  s(t)  and  n(t)  are  equal  and 
given  by 


f  1  „  /  2nWs  _  .  /  ,  2rtW\ 

2W  «  f0r  K  *  T >  S  "  S  <-“o  +  T >  ■ 


V“>  =  p„(“>  -  {  h  > for  ("o  -  s  “  s  (“0  +  . 


0  elsewhere  . 


(A. 3) 
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The  correlation  functions  of  both  s(t)  and  n(t)  are, 
therefore, 


<i>  (t)  = 

ssv 


nn 


*  h 


P  (u)  e^UT  dw 
sx  ' 


sin  jtWt 

1  COS  CJ  T. 

jtWt  o 


(A. 4) 


Also, 


$  a(t) 
ss'  ' 


*  t  \  sin  jcWt 

"nfl(T)  *  — —  sln  V  > 


(A. 5) 


♦h^  (t)  " 


‘»2N2(t)=  *Z1E1  't'  *E2E2  ll‘l  '  itWT  ’ 


^NN  ^  ~  ^  ~  ^  , 

12 


E  E 
12 


(A.7) 


l>  (t)  =  a  $  (t) 

XX  v  '  n  nn  '  ' 


(A.8) 


and 


<t>  (t)  =  <t>  (t)  . 

yy  ss  ' 


(A. 9) 


For  the  regular  correlation  function,  crR  from  Equation  (2.40)  is 
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where  R 
T 
M 

and 

f 

s 


f  /2w  =  ratio  of  sampling  rate  to  Shannon  rate, 
s 

signal  duration, 

total  samples  in  T  seconds 

l/T  -  sampling  rate, 
s 


Equation  (A. 11)  is  plotted  in  Figure  2.6  for  f  /W  =  9  and 
2TW  =  50. 

2 

For  the  quadrature  correlation,  a.  from  Equation  (2.1+6)  is 
given  by 


n 


1 

R 


M/2 

Urn  21 1 

2 

(1  -  — ) 

V  M  ’ 

1  +  2  ) 

R 

m 

m-1 

R 

— 

(A. 12) 


'  2TW1 


where  R  =  f'  /W  and  £  =  l/2T  .  Equation  (A. 12)  is  plotted  in 

s  s  s 

Figure  2.7  for  f  /W  -  9  and  2TW  =  50. 

0  MW 

2 

For  the  clipped  signal  case,  from  Equation  (2.59) »  0.  for 

Q 


quadrature  correlation  is 

WM 
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Ar  'i  {■  +  ([•■■- 

(_  m=l  v 
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R  o 
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(A. 13) 


as  shown  in  Figure  2.9.  The  corresponding  regular  correlation  function 


variance 
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/VWv\ 

2 

gR  =  1 

(1/2TW)  "  R 


(A.lh) 


is  plotted  in  Figure  2.10. 
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